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ABSTRACT 

We present Keck spectroscopic observations and redshifts for a sample of 767 Hersche^SPiRE se- 
lected galaxies (HSGs) at 250, 350, and 500 /im, taken with the Keck I Low Resolution Imaging 
Spectrometer (LRIS) and the Keck II DEep Imaging Multi-Object Spectrograph (DEIMOS). The 
redshift distribution of these Spire sources from the Herschel Multitiered Extragalactic Survey (Her- 
MES) peaks at z = 0.85, with 731 sources at z < 2 and a tail of sources out to z ~ 5. We measure more 
significant disagreement between photometric and spectroscopic redshifts ((Az/(l+z spcc ))=0.29) than 
is seen in non-infrared selected samples, likely due to enhanced star formation rates and dust obscura- 
tion in infrared-selected galaxies. The infrared data are used to directly measure integrated infrared 
luminosities and dust temperatures independent of radio or 24 fim flux densities. By probing the dust 
spectral energy distribution (SED) at its peak, we estimate that the vast majority (72-83%) of z < 2 
_ffersc/ieZ-selected galaxies would drop out of traditional submillimetcr surveys at 0.85-1 mm. We find 
that dust temperature traces infrared luminosity, due in part to the Spire wavelength selection biases, 
and partially from physical effects. As a result, we measure no significant trend in Spire color with 
redshift; if dust temperature were independent of luminosity or redshift, a trend in Spire color would 
be expected. Composite infrared SEDs are constructed as a function of infrared luminosity, showing 
the increase in dust temperature with luminosity, and subtle change in near-infrared and mid-infrared 
spectral properties. Moderate evolution in the far-infrared (FIR) /radio correlation is measured for 
this partially radio-selected sample, with giR oc (1 + z )-0- 30 ± - 02 at z < 2. We estimate the luminosity 
function and implied star- format ion rate density contribution of HSGs at z < 1.6 and find overall 
agreement with work based on 24 /im extrapolations of the LIRG, ULIRG and total infrared contribu- 
tions. This work significantly increased the number of spectroscopically confirmed infrared-luminous 
galaxies at z > and demonstrates the growing importance of dusty starbursts for galaxy evolution 
studies and the build-up of stellar mass throughout cosmic time. 

Subject headings: galaxies: evolution — galaxies: high-redshift — galaxies: infrared — galaxies: star- 
bursts — submillimetre: galaxies 
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1. INTRODUCTION 

Ultraluminous Infrared Galaxies (ULIRGs; Lir > 
10 12 L©) exhibit the most extreme star-formation rates 
in the Unive rse (see an overview in Lonsdale, Farrah fc 
Smith||2006|). At early epochs (z> l), UL1RG activity 



contributes significantly to the build-up of stellar mass 
presumably through intense star-forming burs ts (with 
r < 100 Myr and SFR > SOOM^yr- 1 , e.g. see|Sanders 
fc Mirabel|fl996| |Blain et aZ.~||2002| |Smail et 0^ 2002). 



Since the observed properties of these starbursts are 
short-lived and intense, they are thought t o be triggered 
by th e collision of gas-rich disk galaxies ( |Engel et al. 
2010 1 and serve as a fundamental transition phase to lu- 



minous a£tive galactic nuclei (AGN) or quasars (| Sanders 



et al. ||1988| . Although the merger history of high- z 
ULIRGs has recently come into question, with some evi- 
dence pointing to a substantial (perhaps >50%) fraction 
of ULIRGs building stell ar mass through m i nor mergers 
or passive gas accretion jDaddi et aT]|2010| Elbaz et al. 
|2011| |Rodighiero et al. ||201ip , there is little doubt that 
ULIRGs contribute non-negligibly to the star formation 
history of the Universe and the format ion of massive el- 
liptic al galaxies at the present day (e.g. Kartaltepe et al. 
20101. Unfortunately, much about the infrared star burst 
population 28 is still unknown due to limitations in far- 
infrared (FIR) observations, strong selection biases and 
sample inhomogeneity. 

Galaxies which have been called 'Submillimetre Galax- 
ies' (SMGs), are selected at wavelengths around 1mm, 
particularly in the atmospheric wind ow at 850 um. Such 



'classical SMGs' with S S50 £ 5mJy (Smail et al. 20021 



have put powerful constraints on galaxy evolution theo 
ries and the environments of heav y star formation since 



their initial discovery a decade ago (Smail, Ivison fc Blain 
1997[|Hughes et al. |1998{ |Barger et al. |1999p . Ha 



their selection at 



Dwcvcr, 



tible to strong temperature biasing ( 


Blain et al. 


2004 


Chapman et al. ||2004[ |Case 
al. |2010||Magdis et al. |2010 


y et al. 
1. This le 


2009; |Chapman et 
aves the possibility 



be discovered and characterised. Building a comprehen- 
sive sample of spectroscopically-confirmed submillimetre 
galaxies or extreme starbursts is paramount for deter- 
mining the evolutionary histories of ULIRGs by break- 
ing the Td ust /(1 + z) degeneracy, for carrying out stellar 
population analysis, and for measuring the AGN stage 

and contribution to luminosity. 

The Herschel Space Observatory ( Pilbratt et al. |2010 | 



has identified thousands of galaxies at 70-500 /zm, wave 
lengths previously near-inaccessible from the ground, 
sampling galaxies' emission at the peak of their SEDs 
at z~ 2, when the importance of dusty starbursts in the 
global c ontext of the Universe' s star formation is most 
evident (Chapman et al. 2005[ ). The Spire instrument 



( |Griffin eTai. ||i>0l0| ) will map -350 deg 2 in sky near the 



NSW 1710, Australia 

1 Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with 
important participation from NASA. 

28 In this paper, we use the term 'starburst' to refer to high- 
SFR galaxies (SFR>100 MQyr" 1 ). This differs from the recent 
defi nition of 'starburst' as a combination of SFR and stellar mass 
(e.g. |Noeske et al. [2007] |Rodighiero et al. |2011[ |. 



confusion limit at 250, 350 and 500 /xm as part of the Her- 
schel Multi-tiered Extragalactic Survey (HerMES; Oliver 
et al. in press), covering areas much larger than SCUBA, 
MAMBO, AzTEC or LABOCA. With much larger ar- 
eas, the rarest sources can be uncovered and the dy- 
namic range of sources thereby expands, from nearby lu- 
minous infrared galaxies (LIRGs; > 10 11 L Q ) to distant 
hyper-luminous infrared galaxies (HyLIRGs; >10 13 L Q ) 
and lensed sources. Working towards completeness in 
high redshift starburst samples, and removing the impact 
of selection biases introduced by prior starburst selection 
techniques, is a key long term goal. 

Understanding ULIRG completeness and, in turn, con- 
straining key astrophysical quantities of the luminous 
starburst population, is only possible with a spectro- 
scopic census of a diverse population of FIR-luminous 
galaxies. Redshift identification is a crucial piece of in- 
formation for a high-redshift dusty galaxy, since it allows 
the measurement of its luminosity and star-formation 
rate, and is a prelude to subsequent interferometry (of- 
ten dependent on a known redshift) in order to constrain 
the vast reservoirs of molecular gas which fuel extreme 
starbursts. Such subsequent studies cannot be completed 
using photometric data alone, and to date, Keck LRIS 
and DEIMOS multi-slit spectroscopy is the most efficient 
method for uncovering large samples of gal axy redshifts 
for both normal star-forming z > 1 galaxies ( Cowie fe Hu 



1998) |Cowie, Songaila fc Barger|19991|Cowie et al. |2001| 



Steidel et al. _ 1996[ 1999D and heavily dust-obscured 



ULIRGs dBarger et al. tf5M8\ [19991 |Barger, Cowie~fc" 



Richar ds 2000; |Cowie, Barger fc Kne ib 2002; |Chapman 



et al. i;2005) 



This paper presents the first results from a large 
spectroscopic redshift survey of 1594 Herschel-SPIRE se- 
lected galaxies (HSGs). We measure redshifts for 767 of 
1594 targeted HSGs, describe their bulk infrared prop- 
erties, address their relationship to the now well-studied 
SMGs, and assess their contribution to cosmic star for- 
mation. The results of this spectroscopic survey have 
been split into two papers of which this is the first, 
presenting the details of source selection, completeness, 
spectroscopic confirmations and associated results for 
z < 2 sources. An accompanying paper presents the 
2 < z < 5 sub-sample in m ore detail. Througho ut we 
use a flat ACDM cosmology ( Hinshaw et al. |[2009 ) with 
#0=71 km s" 1 Mpc" 1 and 



^277 



2. HERSCHEL-SELECTEB GALAXY SAMPLE 

The sources observed in this paper were detected by 
the Herschel Space Observatory Spire instrument as 
part of the Herschel Multi-tiered Extragalactic Survey 
(HerMES; Oliver et al. in prep.) . Spire, the Spect ral and 
Photometric Imaging Receiver ( |Griffin et al. 2010 1, is de- 
signed for wide-field mapping at 250, 350, and 500 tun. 
The beamsizes at these respective wavelengths are 18", 
25", and 36" with measured mean point-source confusion 
noise uncertainties of 0250 = 3.8 mJy, 0350 = 4.6 mJy, 
and (7500 =5.2mJy which dominates over instrumen- 
tal noise (these values are for 3cr con f cuts used for de- 
boosted photometric measurements, see t he Spire Ob- 
servers ' Manual and |Griffin et al. 2010 Nguyen et al. 
2010[). We make use of Spire maps as described by 



Levenson et al. (20101 



In this paper, deep ancillary data, particularly radio 



Redshift Survey of Spire sources 
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and 24 /im, are essential for optical spectroscopic sur- 
veying. We observe sources in the Lockman Hole North 
(LHN) whose Spitzer imaging comes from the Spitzer 



dale et al. 



Wide- Area Infrare d Extragalatic (SWIRE) survey Wxms 



2003[ ) and GO MIPS programs (PI Uwen 
and very dee p T.4 GHz mapping from the Very Large 
Array (VLA) (|Owen fc Morrison[[2008l) . LHN has addi- 



tional coverage with Herschel-FACS from HerMES (PI G. 
Magdis). In the Great Observatories Origins Deep Sur- 
vey North (GOODS-N) fiel d, deep 1.4 GHz radio map- 
ping comes from the VLA ( jMorrison et al. |2010[ ) and 
Spitzer coverage of the GOODS-JN center is from FI- 
DEL (Dickinson etal. in prep.) and some of the ex- 
tended area is from Spitzer program ID83 (PI Rieke; 
Shupe, private communication). We also observe sources 
in the ELAIS-N1 (EN1) and extended UKIDSS Ultra 
Deep Field (UDS)/XMM fields, both extragalactic areas 
in the SWIRE survey. The UDS has additional coverage 
from the Spitzer Legacy Program (SpUDS; PI Dunlop). 
Radio coverage of the EN1 field is substantially more 
sparse than in GOODS-N or LHN, wi th the only mapping 
taken with the GM RT at 610 MHz dGar n fc Alexander 



2008[ ) and 325 MHz ( |Sirothia et al |2009p , whose depths 



effectively translate to 100 /zJy RMS at 1.4 GHz assu ming 



a synchrotro n slope of a=0.75, where S oc v a (e.g. Dale 
et al. | |20"07 ). Our UDS observations sit on the edge of a 



new deep VLA radio map with ^7 /iJy RMS (Arumugam 
et al. in pre p.). In the Cosmic Evolution Survey field 



(COSMOS; |Scoville et al 
central 1 de g 2 has a depth o 



20071, radio cove rage in the 
^T0.5 /iJy RMS flSchinnerer] 
et al. |2007[) and Spitzer coverage is described inlanders 



et al. |(200 



(2009 



), fLe Floc'h et al 
MOS is also covered by 
EP program (PI D. Lutz). We additionally observed 
sources in the Extended Chandra Deep Field South or 



and Frayer et al, 
3 ACS as part of the 



ECDF-S region which has dee p radio coverage ( Miller et 
al. ||2008| |Biggs et al. ||2011|), as well as Spitzer-M Lfti 



24 /im Tfom FIDEL and 1RAC ( [Damen et aLpm) . 



We use the photo metric redshift catalogs fr om SWIRE 
(EN1) described by |Rowan-Robinson et al. | ( |2008| , the 
deep LHN catalog in Strazzullo et al. (2010), and the ex- 



tensive COSMOS jnbert et al. ||20l0[ ) a:ul ECDFS (Oar- 



damone et al. ||2010| ) catalogs. In GOODS-N, since we 



survey radio galaxies outside of the central deep region, 
the photometric catalog is limited, so we exclude it from 
photometric redshift analysis. 

2.1. Source Extraction and Photometry 

Spire point source photometry is performed by flux 
extraction at positions of known 24 /zm sources or radio 
1.4 GHz sources. This cross-identification prior source 
extraction (X I D) me thod is described in detail in Rose- 
boom et al. (|2010) wit h some follow-up discussion m 



Roseboom et al. "(2012). The disadvantage of the XID 
technique is that it excludes any sources which are not 
24 /im or radio identified. This is particularly problem- 
atic for potentially high-rcdshift sources which drop out 
of both 24 /im and radio surveys, and might have ambigu- 
ous near-IR counterparts; since 24/iin and radio source 
dropouts are excluded from this sample, and their in- 
fluence is undoubt ably more significant at high-redshifts 
than at z < 2 (e.g. Ma gdis et al. ||20lT] Bethermin et al. 



2012 ), the high-z s ample is treated in a separate paper 



The advantage of cross-identification with 
24/im/IRAC and radio sources is that it can cor- 
rect for confusion boosting in the extracted Spire flux 
densities by estimating the flux contributions from 
nearby sources within one beamsize. It also reduces the 
confusion noise by a factor of ^2 by pushing slightly 
below the nominal confusion limit using the Least 
Absolute Shrinkage and Selection Operator (LASSO) 
method to assign Spire flux densities to an over dense 
list 



prior source 
method details). 



see 



Roseboom et al. 2010 for 



The LASSO algorithm combines 



strengths of model prediction and filter prediction 
source identification (i.e. balancing source priors to the 
Spire map flux distribution with choosing the brightest 
source correspondence between 24 /iin/radio and Spire). 
The algorithm upweights 'rare' sources, therefore radio 
sources are prefered Spire counterparts over 24 /im 
sources. This makes sense giv en the expectation that 



radio sources are FIR-luminou s (Helou, Soifer & Rowan- 
Robinson 1985 Condon 1992 1. The procedure assumes 
that the ancillary data are of adequate depth to identify 
the vast majority of FIR emitting sources, so it is only 
practical in deep legacy survey fields. 

2.2. Completeness of Source Catalog 

The most crucial aspect of a redshift survey is having 
a clear understanding of survey completeness and biases. 
This subsection addresses our survey's completeness in 
identifing strong Spire sources using the XID flux ex- 
traction method. 

Since this is a test of the robustness of the XID tech- 
nique, the results will vary by field (LHN, GOODS- 
N, EN1, COSMOS, CDFS and UDS), based on survey 
depths, and i s a fu nction of Spire flux density. |Rose-| 



boom et al. (20101 demonstrate that the robustness of 
the XID catalog depends greatly on the depth of ancillary 
data available in the field to act as source priors. Beyond 
a fairly standard depth at 24 /tm of S 2 a ~ 150 /iJy (or a 
sky density of ^ 3000 dcg~ 2 ), XID source flux extrac- 
tion will be >95% complete at the 3<r limit of Spire 
(where a includes instrumental and confusion noise, and 
3er roughly corresponds to a ^ 15mJy cut-off). At 
higher source densities XID is more robust, even when 
the source density of priors exceeds the number of Spire- 
bright sources. 

The completeness of the XID process is tested by first 
producing residual maps at the Spire wavelengths using 
the best fit XID solution, and then re-injecting sources 
into these maps using a number count and clustering 
model consistent with the real data. The XID process is 
repeated on these simulated images and the results are 
assessed to determine the number of sources returned at 
> 3er as a function of injected flux density. This process 
is repeated several times for different realisations of the 
Spire maps to build up suitable statistics across a wide 
range of Spire flux densities. The results of this process 
are shown in Figure[T]for the field with the shallowest an- 
cillary data, EN1; deeper fields have completeness curves 
approaching boxcar functions. This completeness curve 



(Casey et ar||2012 henceforth C12) 



Roseboom et 
Roseboom et 



is generat ed using the same analysis from 
al. | fl2010[) with improvements described in 
^pOT2| . An Sspire > 15mJy limit should be >807o 
complete it 500 /tm (>95% for 250 /im and 350 /im), and 
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Fig. 1. — XID catalog completeness curves in the EN1 field 
as a function of flux density at Spire wavelengths 250 /im (blue), 
350 fim (green) and 500 fim (red). The dotted vertical lines mark 
the 3cr noise limits at the respective wavelengths. This complete- 
ness measure is a reflection of the XID catalog's ability to assign 
and identify 24 (im or radio counterparts for all Spire sources at a 
given flux density. 

a 20mJy limit should be >95% complete in all fields. 
Additional sources are identified down to flux densities 
comparable to the confusion limits (5-6 mJy) . This is 
below the nominal 3<r cut-off and is achieved by using 
positional priors of many galaxies thought only to emit 
at Spire wavelengths at levels ~ 0.1-1 mJy; the flux den- 
sity of SPiRE-bright sources is then 'deboosted' using the 
density of source priors which are thought to be Spire- 
faint. 

Note that the ability to match counterparts down 
to 15mJy does not necessarily mean that counterpart 
matching is always correct (an issue which is so far largely 
unconstrained, but is starting to be addressed through 
systematic interferometric work, e.g. Wang et al. [2011 ). 
This question, particularly as it applies to the application 
of the XID extract ion technique, is discussed in |Rose-| 
boom et al. M2010I). 



2.3. Completeness in Radio, 24 fJ-m Samples 

The source density of >3cr SPiRE-significant XID 
sources is constant within a factor of ^2 (1500- 
3000 deg -2 ) despite large variations in ancillary data 
across the different fields. For example, the Lockman 
Hole North — a field with ultra-deep radio coverage and 
deep Spitzer mid-infrared coverage — has a source den- 
sity of ^3000 deg -2 for sources above Si. 4 ~25/xJy or 
S*24 ^150 /iJy, which is the same density measured in 
ELAIS-N1 — a field with only very shallow radio cover- 
age — and in the GOODS-N flanking fields — an area 
not completely covered by Spitzer. Figure [2] shows the 
cumulative surface densities of sources selected at 24 fim 
and radio, then selected to be >3a in at least one of 
the Spire bands. The surface density of 24 /mi Spire- 
significant sources is about a factor of 2 times the sur- 
face density of radio SPiRE-significant sources, assum- 
ing a rough correlation between 24 /jm and 1.4 GHz ra- 
dio flux density of S24 ~ 3.7 Si. 4 (which we measure 
from Figure [3]). Figure [2] also highlights that the sur- 
face density of any one field is dependent on the depth 
of the prior catalog, although across all fields this is 
consistent within a factor of about 2. With compara- 
ble surface densities ~ 1500-3000 deg -2 , one might then 
ask whether this means that 24 /im identified sources 
in EN1/COSMOS/CDFS/UDS and the radio-identified 
(>40/xJy) sources in LHN/GOODS-N/COSMOS/CDFS 



are drawn from the same population of IR galaxies. 

Figure [3] investigates the relation between 24 fim and 
1.4 GHz flux density for the LHN subsample which con- 
tains the deepest ancillary data. Our Spire targets are 
shown against normal field galaxies, both carving out 
a similar parameter space, indicating no clear Spire- 
bright/faint bias. Using our LHN sample as a guide, 
we find that 70±20% of Si. 4 > 40/xJy sources have 
S 2 4 > 150/iJy (see Fig EL suggesting that most sub-mJy 
radio sources are sub-mJy 24/im sources. This highlights 
that: (a) the incidence of SPiRE-detection is higher in 
faint radio sources than in 24/xm sources; and (b) that 
nearly all radio sources (^70%) are 24 (im detected and 
a large fraction of 24 fim sources are radio detected (from 
Fig 4h. Since our sample is rest-frame FI R selected and 
is likely to obey the FIR/ radio correlation ( |Helou, Soifer| 
fc Rowan-Robinson|1985 ), our sources are more likely to 



be radio-detected than random field galaxies. Also, given 
the beamsize of Herschel to be 18-36", the probability 
of random coincidence with a radio galaxy is >8 times 
lower than with a 24/im galaxy, thus the probability of 
correct counterpart identification is >8 times higher in 
radio galaxies than it is in 24/im galaxies (although that 
probability itself cannot be constrained without interfer- 
ometric infrared observations). 

The typical range of surface densities of Spire galax- 
ies, ^1500-3000 deg -2 , translates to an expected number 
of SPiRE-detected sources per slit-mask of Ndeimos = 
33^g 5 sources per mask and Nlris = 17+* sources per 
mask across all fields (derived from mean and standard 
deviations of source densities between the six survey 
fields). This agrees with our actual spectroscopic sam- 
pling per slit-mask; we average 33±7 significant Spire- 
sources per DEIMOS slit-mask and 20±6 sources per 
LRIS slit-mask. 

It is important to note that this analysis only tests 
the completeness and source density for galaxies already 
detected at mid-infrared or radio wavelengths. Some 
sources will be excluded from the XID catalog since they 
will not be detected in the ancillary data. Important 
examples are very high redshift sources which are radio- 
faint and 24 /im-faint and cannot be identified due to a 
lack of a multi- wavelength counterpart (s). These high- 
redshift sources will have a profound effect on the derived 
IR luminosity function at z > 3; therefore, this is a signif- 
icant limitation of this HSG sample. However, at lower 
redshift we suspect that, at the least, the XID catalog 
is 80% complete out to z = 2 (again, see more det ails on 
estimated complet eness as a function of redshift in |Rose- 



boom et al. |2010 ). The completeness of the XID catalog 



with respect to the low redshift sources, z ^ 2, has also 



been shown to be «95% ( |Magdis et al. ||2011[ ) . 

2.4. Spectroscopic Target Characteristics 

Slit-masks for t he Keck Low Resolu tion Imaging Spec- 
trometer (LRIS; Oke et aLj [l995) and DEep Imag- 



n 



i ng Mu lti-O.jbect Spectrograph (DEIMOS; Faber et al. 

prioritization scheme 



2003 ) were populated by a 
whereby sources were graded in priority from to 1000 
by their Spire photometry, radio detection or non- 
detection, and 24 /im flux density. No optical magnitude 
selection or prioritization was made. Note that there was 
no filtering of the sample to remove low re dshift targets 
or quasars (as was often done for SMGs, e.g. Chapman et 
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Fig. 2. — Cumulative surface density of 24 and radio sources with flux density. In other words, as a function of radio or 24[ixa 



flux density, wc plot the surface density of sources with flux densities above that given flux density. Dashed lines represent the parent 
population of all 24 am sources (blue) and radio sources (red) in the whole held. Each panel is labeled at top with a 24/im flux density 
scale and at bottom with a 1.4 GHz flux density scale (the scales arc omitted it the data arc insufficient/do not exist). Solid lines represent 
HTe"^ouree"'clCTTsTt^ToTl5oTirce^^ 

significant in all three Spire bands. Vertical dotted lines mark where 24 ara or radio catalogs become incomplete in flux density. 



al. |2005[ ). The i-band magnitudes of our spectroscopic 



targets are plotted against the brightest Spire flux in 
Figure [5] Since our spectroscopic sampling is exclusively 
driven by Spire detectability and not by optical magni- 
tude, there is no clear relationship perceived in Figure [5j 
When there was no «-band or z-band counterpart to cen- 
ter our slit on, we used the IRAC 3.6 /im position or radio 
VLA position, which are both good to ^0.7" given their 
relatively small bcamsizes. 

Sources detected at >3tr in at least one out of the three 
Spire bands were given a priority = 300, and sources 
detected in all three Spire bands at >3<r were given a 
priority = 500. Mask centers and orientations were cho- 
sen based on the positions of rare, 'red' 500 /xm-peaking 
sources (S^so < S350 < •S'soo), t hought to be the highest- 
redshift SPiRE-bright galaxies ( |Cox et a7~||2011 ). These 
red sources were given very high priority, >800, and if 
their multi-wavelength properties were consistent with a 
high-z source, e.g. zab >22 and S 2 a <500 /iJy, prior- 
ity was graded even higher, at 1000. The prioritization 
scheme is linear, such that a source with priority 1000 
will be assigned a slit in favor of two p=300 sources, or 
one p=300 source plus one p =500 source, however, a 
p=800 source would be passed up in lieu of two p=500 
sources. While this scheme could accidentally remove 



very high priority sources from our masks, we adjusted 
mask position angles manually to ensure optimal spatial 
sampling, and were only minimally affected by source 
overlap ( 5% of slits were conflicted). 

The density of our sources, both low and high priority, 
is high enough and comparatively uniform over the LHN, 
GOODS-N, ENl, UDS, CDFS and COSMOS fields such 
that our mask coverages constitute a random sampling, 
and completeness estimates are performed with respect 
to the sky area probed by the masks alone (which to- 
tal ^0.93 deg 2 ). Note that the centering of our masks 
around high-priority 'red' sources might give a high- 
redshift bias to our sample, since we set out to find some 
of the rarest, high-z HyLIRGs. What we find (to be 
discussed later in the paper and more in C12) is that 
the redshift distribution of 'red' targets is not strongly 
biased towards high-z. There are actually more high-z 
sources which are not red than are, leading us to believe 
that the sky sampling is essentially random, despite our 
efforts to detect more high-z galaxies. The density of 
>3cr Spire sources is also low enough so that our radio, 
24 /im or Spire color prioritization does not introduce 
statistically significant selection biases into the sample; 
in other words, nearly all (>90%) the >3cr Spire sources 
within areas covered by slit-masks were spectroscopically 



G 



C. M. Casey et al. 



10000 



1000 



100 




24 



1000 



Fig. 3. — The distribution of LHN targets in 1.4 GHz radio flux 
density (Si. 4) versus 24 /im flux density (S l 24). The various flux 
limits in the three different fields are marked. Black points repre- 
sent field galaxies while red points denote SPlRE-selected galaxies. 
Note that Spire detection is largely uncorrelated with 24 (im or ra- 
dio flux density. However, faint 24 fim sources tend also to be faint 
in the radio (roughly corresponding to S24 ~ 3.7 X Si. 4 as mea- 
sured here), leading us to conclude that the majority of faint 24 (im 
EN1 sources are drawn from the sa me population as the G OODS-N 
radio sources (also see the work of |Magdis et al. |2011| . 
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Fig. 4. — The fraction of radio sources (red) or 24 fim sources 
(blue) which are SPlRE-detected as a function of radio or 24 fim 
flux density. The 'x's represent the fraction of 2Afim Spire sources 
which are also radio-detected (blue) and radio Spire sources which 
are also 24/xm-detected. Here we use S24 ~ 3.7 X Si. 4, which is the 
rough relation we observe for field galaxies in Figure [3] this scaling 
allows us to judge the relative depths of the catalogs. 'Depth A' 
corresponds roughly to Si. 4 ~40/xJy and S24 ~150/xJy, 'Depth 
B' corresponds roughly to Si. 4 ~25 fijy and S24 ~100^tJy, and 
'Depth C' corresponds roughly to Si. 4 ~15/^Jy and S24 ~70/^Jy. 
'Depth A' corresponds to the EN1 and UDS 24 fim limit, 'Depth B' 
corresponds to LHN 24/im and COSMOS/GOODS-N/CDFS radio 
limits, while 'Depth C corresponds to the LHN radio limit and 
GOODS-N(center)/CDFS 24 fim limit. 

observed. 

We assigned slits for 100% of the very high priority 
'red' targets, 95% of the high priority targets (>500), and 
90% of the lower priority targets (300< p <500), filling 
any free space with additional objects not discussed in 
this paper (radio sources or 24 fim galaxies with insignif- 
icant Spire fluxes). Any high-priority Spire sources not 
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Fig. 5. — Optical i-band magnitude (AB) against Spire flux 
density (the brightest flux density is taken from Spire 250 fim, 
350 fim or 500 fim) for sources in EN1 (similar results for other 
fields). Black points are all significant Spire sources used to fill our 
slitmasks. The selection of these Spire sources is not a function of 
optical magnitude. 

observed were only excluded on the basis of the mask con- 
figuration, where observing them would bump another 
high-priority source off a slit. However, the exclusion of 
high-priority Spire sources is rare enough that having a 
few missing from our sample does not impact our com- 
pleteness, particularly since we find no significant bias 
towards high-redshift spectroscopic identifications in the 
high-priority target sub-sample (see C12 for more details 
on the z > 2 sub-sample). 

Our targets' distribution in radio flux density and in 
24 fim flux density and Spire signal-to-noise is shown in 
Figure [6j Out of 1594 targets, 633 are 1.4 GHz identified, 
1319 are 24 fim identified, and 588 are both radio and 
24 fim identified. The sources which are selected in both 
radio and 24 fim maps come primarily from the LHN, 
COSMOS or CDF-S fields, since EN1 and UDS lacks 
deep radio coverage and many sources in the flanking 
fields of GOODS-N lack deep Spitzer 24 fim coverage. 

Throughout this paper, we refer to our sample as Her- 
schel SPlRE-selected galaxies (HSGs). We prefer not to 
call them submillimetre galaxies (SMGs) for sake of con- 
fusion with the historical definition of SMC selected at 
850 fim— 1.4 mm. We also choose not to use the term 
ULIRG for our sources since it places the special qual- 
ification of a luminosity cut on the sample (10 12 L Q < 
L < 10 13 L Q ). For the rest of the paper, any reference 
to SMGs refers to a population selected at 850 fim with 
5*850 >5mJy. We use this strict definition of the SMG 
population in order to draw comparisons between 'classic 
SMGs' and HSGs. 

3. SPECTRAL OBSERVATIONS AND REDSHIFT 
IDENTIFICATION 

Optical spectroscopic observations were carried out at 
the W.M. Keck Observatory using LRIS on Keck I and 
DEIMOS on Keck II in 2011 and 2012. 

LRIS observations were carried out in adequate condi- 
tions on 2011-Feb-06 with ~1 " seeing and cloud cover, on 
2012-Jan-26 and 2012-Jan-27 with -0.8-1.2 " seeing with 
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Fig. 6. — Top: Distribution of 24 /im and 1.4 GHz radio flux 
density of our targets. Of our sample, 1319 were primarily identi- 
fied through 24 fim emission and 633 of were radio selected. The 
gray distributions reflect the sources with confirmed spectroscopic 
redshifts (this number exceeds 767 since some sources are identi- 
fied as both 24 fim and radio sources). The spectroscopic success 
rates for both samples is ~55-60% and does not correlate with 
flux density. Dashed vertical lines represent average flux density 
limits of the radio data (~30 fiJy) and the 24 (im data (~150 fiJy) 
in our fields. Bottom: Each target in our survey must be de- 
tected at >3cr in at least one of the three Spire bands at 250, 
350, or 500 fim. Here we show the 250 fim signal-to-noise ratio 
(SNR) against 350 /im SNR. The 500 fan SNR scales with the size 
of the gray circles; if a target is under 3<r at 500 /im, it has no gray 
circle. Spectroscopically confirmed sources are then marked with 
red crosses, showing no obvious correspondence between FIR flux 
density and spectroscopic success rates. 

minor to no cloud cover, and on 2012-Feb-27 with —0.5- 
1.0" seeing and no cloud cover. We used the 400/3400 
grism for maximum wavelength coverage in the blue. 
The 2011 observations used the 600/7500 grating in the 
red with a central 'multi-slit' wavelength of 6500 A; the 
2012- Jan observations used the 400/8500 grating in the 
red with a central wavelength of 8400 A, and the 2012- 
Feb observations used a central wavelength of 8000 A. 
All observations used clear filters for both red and blue 
arms and the 560 nm dichroic. These setups give a 1.09A 
dispersion in the blue (e.g. shortward of 5600 A) and a 
0.80 A dispersion in the red. Wavelength coverage for 
each source varied by its position on the 5.5' x 7.8' slit- 
mask, the mean wavelength coverage ranged from 2500- 
8200 A in 2011 and from 2500A-1 //m in 2012 but varied 
up to —1300 A for sources on each mask. Due to the 
dichroic, some sources near the edge of the slit-mask have 
gaps in wavelength coverage -800 A wide in the vicinity 



of 5600 A. 

We observed a total of 25 LRIS multi-slit masks, 
13 of which were observed in near-photometric condi- 
tions. LRIS data reduction, including bias subtraction, 
flat fielding, wavelength calibration, and sky subtraction 
were all performed using custom-built IDL routines. Of 
664 LRIS targets with >3cr Spire detections, 268 were 
spectroscopically identified. 198 of the identified sources 
(74%) were confirmed in near-ideal conditions; we add 
asterisks to the NAME column of Table [2] to distinguish 
these sources from sources confirmed in poorer weather. 

DEIMOS observations were carried out in good condi- 
tions on 2011-May-28 and 2011-May-29 with -0.6-0.7" 
seeing, on 2011-Nov-28 in average to cloudy conditions 
with -1.0-1.3" seeing, and on 2012-Feb-16 and 2012- 
Feb-17 in very cloudy conditions with 1.0-3.0" seeing; 
we used the 600 lines mm -1 grating with a 7200 A blaze 
angle (resulting in dispersion of 0.65 A) and the GG455 
filter to block out higher-order light. Wavelength cover- 
age varied with source position on the 5 ' x 16.7 ' slit-mask 
from 4400-9200 A to 5200-9900 A and averaged to 4850- 
9550 A. 

Sixteen of 29 DEIMOS multi-slit masks were observed 
in near-photometric conditions, with integration times 
- 2700-4800 s. We used the DEEP2 DEIMOS data re- 
duction pipeline to reduce these data 29 . Of 930 DEIMOS 
HSG targets, 499 were spectroscopically confirmed. 324 
of those 499 (65%) were confirmed in near-ideal condi- 
tions and are also marked with asterisks in Table [21 

Redshift identification was carried out through the 
identification of multiple spectral signatures, primarily 
with the [Oil] doublet, [OIII], H/3, Ha, [NII],CaH&K 
absorption and the Balmer break, H7, Lya , and the Ly- 
man break (given in order of decreasing occurance in the 
sample). A minority of sources were identified by CIII], 
CIV and Hell emission (all sources identified via fea- 
tures in the rest-frame ultraviolet are discussed in C12). 
Full spectroscopic details of all confirmed z < 2 HSGs 
are given in Table [2] (available online in full). The relia- 
bility of redshift identifications is also given in Table [2j 
rated on a scale of g 2 =l-5, where 5 is best. Sources with 
multiple feature identifications have q z > 3, and sources 
with single line identifications have q z — 1-2. When visu- 
ally inspecting the spectra for our targets, redshifts were 
graded on a wider scale, with some sources having q z =0 
(poor) or -1 (non-existent). These sources' potential red- 
shifts are not reported in this paper due to unreliability. 

There are more HSGs with unconfirmed redshifts (826) 
than there are with redshifts (767). The vast majority 
of the unconfirmed sample is unconfirmed due to poor 
weather including cloud-cover and poor seeing (54% of 
sources were observed in sub-optimal conditions). The 
average spectroscopic yield during good weather was 
60±20% while the poor weather yeild was only 15±10%. 
The remaining unconfirmed sources either very red faint 
continuum without identifiable spectral features (20%) 
or no continuum whatsoever (80%). These sources are 
potentially misaligned on the slit (which is based on an 
IRAC and i— or z— band image when available) or are 
too optically obscured to be detected in bright emission 

29 The analysis pipeline used to reduce the DEIMOS data was 
developed at UC Berkeley by Michael Cooper with support from 
NSF grant AST-0071048 
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Fig. 7. — TOP: Distribution in i-band magnitudes of spectro- 
scopically confirmed (hashed area) and unconfirmed sources (filled 
gray area). There is clearly a bias in sources with spectroscopic 
confirmations being brighter at i-band magnitude, also much more 
likely to sit at lower redshifts. However, note that for sources not 
included on this plot — those undetected in i-band — the spectro- 
scopic success rate is 39%, lower than the same fraction for sources 
with i-band detections, 54%, although not as low as one might 
expect for very optically faint sources. This is likely caused by 
many sources being confirmed through bright emission lines which 
can contribute minimally to broad-band photometry. BOTTOM: 
The i-band magnitude against redshift for sources with confirmed 
redshifts. At the bright end, i < 22, sources show a clear trend 
with redshift, while many sources (marked as upper limits, accord- 
ing to the imaging depths in each field) have no available i-band 
photometry. Clearly i-band magnitude strongly impacts the like- 
lihood of measuring a spectroscopic redshift, a concept which is 
explored whe n measuring the spectroscopic completeness, shown 
in Figure [12] 

lines with a ~l-2 hr integration. They could also be 
at z 3, thus intrinsically much more difficult to detect 
optically. 

Since we have very limited constraints on the spectro- 
scopic failures, we emphasise that our sample is spectro- 
scopically incomplete. We use photometric redshi fts a nd 
inferences on the redshift distribution itself in § |3.2| to 
estimate the spectroscopic completeness as a function of 
redshift. 

3.1. Biases in Spectroscopic Confirmations 

Here we quantify the biases of our spectroscopic obser- 
vations by analyzing sources which failed to yield spec- 
troscopic identifications. This is arguably the most dif- 
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Fig. 8. — Relationship between i-band magnitude and Spire 
color, [250] - [500]=— log (S25o/Ssoo)- Spire color becomes redder 
with fainter magnitudes, which is expected given the redshift-i- 
band magnitude relationship seen in Figure Vf\ and an assumption 
of constant dust temperature from low to high redshifts. The lo- 
cus of all >3tr Spire sources in the upper panel consists of ~5000 
sources in the COSMOS field (black dotted contours and black 
points, the mean Spire color for i magnitude bins). Red dot- 
ted contours and points represent the distribution of spectroscopic 
targets (~1600 sources across all fields) and blue represent sources 
with spectroscopic confirmations, clearly skewed towards brighter 
i-band magnitudes. In the bottom panel we compute the mean 
photometric and spectroscopic redshift for the same colored sam- 
ples from the plot above, and see agreement down to iAB ~ 25.5 
(vertical dot-dashed line). At top, the gray band represent the ex- 
pected Spire colors at the mean photometric redshift per i— band 
magnitude bin for dust temperatures ranging from 30 K to 50 K. 
Within uncertainty, the full sample follows this expectation, while 
the spectroscopic sample deviates toward redder colors at brighter 
i-band magnitudes. 

ficult completeness to quantify, since it requires some 
knowledge of the redshift distribution of sources which 
are: (a) the most optically obscured, (b) have feature- 
less continuua, or (c) have observed-frame emission lines 
outside of the wavelength range of our observations (this 



redshift range is 1.6 



3.2 for DEIMOS observations 



and 1.6 ^ z ,$ 1.7 for LRIS observations). Figure[7|shows 
the optical i-band magnitude distributions for sources 
spectroscopically confirmed and unconfirmed, and i-band 
magnitude against redshift. This makes it clear that op- 
tical i-band magnitude need not be very bright for a spec- 
troscopic identification based on emission lines, and that 
optical magnitude itself does not constrain redshift, given 
the number of sources at low- 2 which are undetected in 
the i-band. However, it is clear that there is an overall 
trend with redshift at brighter magnitudes i < 22, with 
very few sources at these magnitudes at z > 0.5. In that 
sense, it is also clear from Figure [7] that optically bright 
sources are more likely to be spectroscopically confirmed 
than those which are faint; the mean i-band magnitude of 
spectroscopically confirmed sources is iAB = 22.1, while 
for unconfirmed source it is iAB = 23.8. 

Prior studies of optical photometric properties of in- 
frared starburst galaxies find poor coorespondence be- 
tween optical magnitude, IR luminosity and redshift (e.g. 



Chapman et al. 2004 2005). This makes it quite diffi 



cult to estimate the redshift distribution and IR hum 
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Fig. 9. — The spectroscopic redshift distribution of HSGs. 
Since the majority of our 767 spectroscopically confirmed sources 
were observed with DEIMOS (66%), the redshift range 2=1.6-3.2 
(hashed area) is sparsely populated (only two sources within this 
range are DEIMOS identified, both bright quasars). Distribution 
of photometric redshifts for the sample of spectroscopically tar- 
geted HSGs is shown as the dotted line, consistent at z ;$ 1.5 but 
suggesting a significantly decreased spectroscopic yield at z >1.5. 

nosities of the optically faint targets for which we fail 
to measure redshifts. Without redshift information, we 
must rely on color information to infer if the optically 
fainter targets have a significantly distinct redshift dis- 
tribution (this also assumes SPiRE-color varies with red- 
shift and dust temperature is roughly fixed). Figure [8] 
shows how Spire color relates to i-band magnitude for 
the entire parent galaxy sample of >3er Spire sources, 
spectroscopic targets, and also spectroscopic confirma- 
tions. This indicates that Spire colors should become 
redder with fainter i-band magnitude, assuming a con- 
stant temperature and redshift range. The mean red- 
shifts — photometric and spectroscopic — per i-band bin 
are self-consistent, indicating no strong bias in the spec- 
troscopic redshift distribution with i magnitude. Note 
however that the analysis as presented in Figure [8] ex- 
cludes sources (a) without photometric redshifts, and (b) 
without i-band counterparts. The variation of Spire 
color with redshift, along with the perceived bias of 
our spectroscopic sample towards redder colors at bright 
magnitudes, is addressed further in § 3.5 

Of the spectroscopically confirmed subsam- 
ple, 65±15% are 250 /im-peaking (meaning 
■5*250 > > 5500), 23±14% are 350 /im-peaking 

(S350 > S250 and #350 > £500), and 13±9% are 
500 ^m-peaking (S250 < S350 < •S'soo)- These values 
are comparable to the spectroscopic targets: 72±12% 
were 250/mi-peaking, 17±9% were 350/mi-peaking and 
11±6% were 500/zm-peaking. This leads us to conclude 
that successful spectroscopic identification does not have 
a bias with respect to FIR SED shape. 

3.2. Redshift Distribution 

The redshift distribution of the 767 spectroscopically 
confirmed HSGs in this survey is shown in Figure [9j Of 
767, the vast majority, 731, are at z < 2. The distri- 
bution peaks at z = 0.85 with a tail of sources extend- 



ing out to higher redshifts, discussed fully in a separate 
accompanying paper, C12. Also plotted is the distribu- 
tion in photometric redshifts for all spectroscopic targets. 
The photometric redshift distribution peaks at the same 
epoch, but with a much higher fraction of sources at 
z ^ 2. The deficit in spectroscopic redshifts at z > 2 is 
caused by: (a) the DEIMOS 'redshift desert' (DEIMOS 
observations comprise 66% of our sample and no strong 
emission lines are visible within DEIMOS wavelength 
coverage from 1.6 < z < 3.2); and (b) by a decreas- 
ing spectroscopic completeness due to enhanced obscu- 
ration in the rest-frame ultraviolet relative to rest-frame 
optical. The latter point relates to the very obscured 
nature of infrared-luminous galaxies: the increased pres- 
ence of dust implies more significant extinction in the 
UV and optical, with more substantial effects at bluer 
wavelengths. 

Note that predicted redshift distributi ons for Herschel 
s ources have been studied in detail b y |Amblard et al.\ 
[] ( |2010[ ) and |Bethermin et al. | ( |2011[ ). However, these 
are based on different SPIRE selection methods and are 
thus not directly comparable to the HSG sample dis- 
cussed herein. For instance, the Amblard et al. distribu- 
tion of 350 /mi-selected sources peaks at z ~ 2.2, as seen 
in the upper panel of Figure 10 We extract HSGs from 
our sample which would satisfy their detection criteria 
(S 350 > 35 mJy, with 250 /an and 500 /im SNR > 3) for 
comparison and find a statistically distinct distribution 
from Amblard et al. 

We also compare data with the model predictions of 



the ex p ected Spire dist r ibutio ns fr om |Lagache et al. 



fl2004| ), |Negrello et al. 

data 



( 2007 ) and |Bethermin et al. 
(T5oth 



(1201 ip . While our data (both of the limited &350 > 
35mJy sample and the distribution of all photometric 
redshifts) are inconsistent with most models (particu- 
larly at the z > 2 end), our results are most consistent 
(tested via a Kolmogo rov-Smirnov statistic) with the pre 



dieted distribution of Lagache et al. (2004); the number 



of sources we observ e at z < 1 is notably d ifferent than 



the predictions from Negrello et al. 
long- wavelength ph otometric- redshift 
|blard et al. Wm^ . 

In the mid dle panel of Figu rc|10|wc compare our results 



(120071 ) or fro m the 
Sasea work of |Am-| 



to those of Chapin et al. (2011) for BL AST-detected 



galaxi es in ECDF-S, a combination of photom etric red 



shifts ( jDunlop et a l. 2010; Iviso n et al. 1 2010 ) and spec- 
troscopic redshifts ( |0asey et al. |2011[ ). Their selection is 
based on three different flux cuts at 250, 350, and 500 ^m, 
and when applying the same selection to our sample, we 
find an agreement between the median redshifts, both 
of which have (z) rj 1 but with differences at both low 
and high-redshift. There is an additional peak in the 
Chapin et al. sample at z < 0.5 which is not observed in 
our sample. This could be due to cosmic variance, since 
the BLAST ECDF-S sample probes a small volume in a 
single field, while our data sample multiple deep fields 
over a larger sampling volume. Given the DEIMOS red- 
shift desert, comparison cannot be drawn fairly beyond 
z~1.6. 

Worth noting is the contrast of these redshift distribu- 
tions with the 850 um-selected SMG redsh ift distribution 



(Chapman et al. 



2005 Yun et al. 



at z 



2012J, which peaks 
see bottom panel of Figure 10). The peak 



in 850 /im-selected galaxies occurs at earlier epochs due 



10 



C. M. Casey et al. 



20 - 



15 



< 

CD 
O 
I— 

cn 
x 




< 

cd 
o 
i— 
cn 

x 



60 - 



40 - 



20 - 



< 

CD 
O 
I — 
CO 

X 



5 Ja0 >35mJy, Our Doto 

Our Dolo (phot-z) ^— ^— 

Amblord+10 Doto 

Logoche+04 Models 

Negrello+07 Models 

Bethermin+ 1 1 Models — — — 



S 250 >40mjy OR Sj5|,>30mjy or S 50l) >20mjy, Our Doto 
Chopin+11 (BLAST) 




All HSGs: Our Doto (spec-z) 
All HSGs: Our Doto (phot-z) 
Chopmon+05 850^m SMGs 
Aretxogo+07 850ju,m SMGs 
Yun+12 1.1mm SMGs 




Fig. 10. — Comparison of our redshift distributions to those in 
the literature; at top, we compare to the Amb lard et al. ] | |2010| l 
distribution (green) for sources selected at 35U /im above 35mjy 
(with >3rx significant detections at both 250 (im and 500 (im). We 
also compare the full distribution of photometric rcdshift s of our 
targets (thick, solid gray) to t he model di stributions from|Lagache| 
\et al. | | |2004) , |Negrello et al. | ( |2007| >, and |Bethermin et a/~ffi01l| 
(dot-dashed red, dashed magenta, and long-dashed blue respec- 
tively). Neither observed distribution (spectroscopic and photo- 
metric) shows good agreement with the models, however the La- 
gache et al. model (which peaks at * ~ 1) agrees best with o ur data. 
The middle p anel compares to the redshift distribution of |Chapin] 
\et al. \$20Uj of for BLAST-bright galaxies in ECDF-S. Given the 
limited number of sources in the BLAST analysis (and that they are 
all identified in ECDF-S, a small volume) we conclude that these 
two distributions are consistent within uncertainty. At bottom we 
compare the HSG redshift distrubiton to those of 850/im-l.lmm 
selected SMGs, which peak at z ss 2.2—2.6 due to the longer wave- 
length selection. The HSG distribution is taken from Figure [9] and 
scaled down 10 times in total number for direct comparison. 



to the selection wavelength: 850 /im is more sensitive to 
lower luminosity/colder sources at higher redshifts and 
less sensitive to warmer sources at low redshift. Spire 
selection is not constant in luminosity out to very high 
redshifts like ~lmm selection is, but it has the benefit 
of probing the peak of the infrared SED at z ~ 1-2, thus 
is not likely to miss or select against sources at a given 
luminosity due to their warmer SED shape. 



3.3. Comparison to photometric redshifts 

Photometric rcdshifts are far easier to obtain on larger 
samples of galaxies than spectroscopic redshifts, and 
whenever the latter become available, it is important to 
test the reliability of the former. Many large statistical 
studies are now motivated exclusively by use of large cat- 
alogs of photometric redshifts; here we explore possible 
underlying biases which might persist in HSG photomet- 
ric redshifts and what might cause them. Note that all 
photometric redshifts used in this paper are calculated 
from ultraviolet through near-infrared photometry and 
exclude any long-wavelength data >8/im (the catalogs 
used in this analysis are described at the end of § [2]) . 

Spectroscopic and photometric redshifts are plotted 
against one another in Figure |11| From this we mea- 
sure that HSG photometric redshifts have a mean rms 
scatter of cr=Az/(l + z sp ec)=0.29, characteristically a 
factor 3-4 times worse th an the photometric redshifts 

measure 



for regular field galaxies 
Az/(l + z spcc )=0.07for ga 



2010 



( Ilbert et al. 

axies with spectral confirma- 
tions in COSMOS). This might be suprising given that 
the majority of sources in this comparison come from 
fields with deep, multi-band photometry (e.g. COSMOS, 
LHN). Why are the photometric redshifts of infrared- 
bright galaxies substantially worse than those of most 
field galaxies? 

One might think that the disagreement originates from 
source blending or mismatching due to the large beam- 
size of infrared observations. However, we do not expect 
blending or mismatching to occur since the photometric 
redshifts taken from the source catalogs are matched to 
our spectroscopic targets' position within 1" (i.e. even 
if a spire source's counterpart is mistakenly identified, 
we would still expect the false counterpart's photometric 
redshift to agree with its spectroscopic redshift). Since 
mismatched counterparts is not likely to be the source of 
the photometric redshift scatter, the scatter is likely due 
to some other property of the infrared-selected sample. 

If the Zphot — z spec disagreement is not caused by mis- 
matched counterparts, then what intrinsic physical pro- 
cesses could explain less reliable photometric redshifts 
in infrared-selected galaxies? Direct detection at far- 
infrared wavelengths implies two things about a galaxy: 
(i) it has a significant dust reservoir which has ab- 
sorbed more energetic light and re-radiated it in the 
far-infrared, and (ii) it has a high star- format ion rate 
(since FIR luminosity scales directly to star formation 

1- 



rate: 



Kennicutt 



1998[ ). The first point highlights that 
the galaxy's emission at rest-frame ultraviolet and op- 
tical wavelengths-emitted by young, hot stars-is being 
obscured and scattered by dust; as a result, infrared- 
selected galaxies are optically fainter than their less- 
dusty counterparts of similar redshifts and stellar masses. 
The second point, that infrared-galaxies have intrinsi- 
cally higher star formation rates than most 'normal' 
field galaxies follows from the implied infrared luminosi- 
ties of the infrared-selected samples. The vast majority 
of field galaxies have star formation rates <lOM0yr _1 
(e.g. the vast array of 200K galaxies in COSMOS which 
are fit with st ellar population templates from |Bruzual fc| 
Charlot|2003 ). In contrast, most infrared-selected galax- 
ies satisfy Ltr > 10 11 L Q , which implies star formation 
rates of > 17M0yr _1 or Ltr > 10 12 L Q which implies 
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Fig. 11. — TOP: Redshift distributions for HSGs shown on 
a logarithmic scale: the full spectroscopic sample (solid black); 
spectroscopically confirmed sources with photometric redshifts 
(zspec I z phot > in other words "z spoc given z p hot>" dashed red); 
spectroscopically confirmed sources without photometric redshifts 



(Zspcc \ -Zphot! in other words "z s 



take ■ 



away z phot ,' 



dotted 



magenta); photometric redshifts (solid blue); and photometric red- 
shifts for spectroscopically confirmed sources (zphotl^spec, in other 
words "Zphot given z BpBC " dashed green). The DEIMOS redshift 
desert is shown as a hashed area. BOTTOM: Spectroscopic ver- 
sus photometric redshift for successfully identified targets which 
have photometric redshifts (43 1 sources out of 767). T he photo- 
metric redshi fts were taken from|Strazzullo et al. |l|2010} (Lockman 
Hole North), JRowan- Robinson et at. | J!jUU8t (ELAIF?N1 UPS) , 
|Hbert et al U'20W\ (COSMOS) and |Cardamone et al. \ \2010\ 
(EUDFbJ with associated uncertainties on z p hot shown. Trie mean 
RMS scatter at is Az/(1 + z spcc )=0.29, notably worse than the 
expected photometric redshift accuracy of photometric redshifts of 
'normal' starburst galaxies. 

> 17OM0yr~ 1 . Higher star formation rates translate 
to brighter rest-optical emission lines and higher line-to- 
continuum ratios. 

How might dust obscuration and enhanced line-to- 
continuum ratios impact the reliability of photometric 
redshifts for HSGs? The first effect is straightforward: 
the fainter a galaxy is in the optical, the more difficult it 
is to put a reliable constraint on its photometric red- 



shift. Since HSGs are optically fainter than 'normal' 
galaxies which harbor less dust, this is one reason HSG 
photometric redshifts are less reliable. The second effect 
comes from the 'contamination' of bright emission lines 
of the optical broad band filters used to compute the 
sources' photometric redshifts. Although emission lines 
can contaminate the broad band magnitudes for both 
normal and dusty galaxies, dusty galaxies have compa- 
rably fainter continuum (or rather, a more significant 
contribution from emission lines). Also, despite some ef- 
forts to account for emission lines in photometric redshift 
code algorithms (by including them in the templates used 
for the fitting), no ne of the existing optical templates in 
the literature (e.g. Calzetti, Kinney k, Storchi-Bergmann 
|1994| [Bruzual fc Charlot||2003p include models with ex- 
tremely high star formation rates and dust content as 
exists in HSGs, hence the poorer z prio t — z svcc agreement. 

The influence of dust obscuration on photometric red- 
shift estimates is not completely straightforward in that 
dust is not necessarily expected to extinct optical flux 
uniformly across all wavelength regimes. Sources with 
high phot-zs and low spec-zs can be explained by dif- 
ferential blue-to-red obscuration whereby the galaxies 
dropout in several blue bands due to stronger absorp- 
tion of higher energy photons. The handful of sources 
with low photometric redshifts and high spectroscopic 
redshifts (z > 2) are described in detail in C12; these 
are thought to disagree due to differential obscuration 
of resonant e mission-line photons and continuum (e.g. 
Neufeld|1991 1 
fi> 



whereby emission lines are not extincted 
as significantly as stellar emission continuum. This is 



also highly dependent on viewing geometry, but the ob- 
servation holds that photometric redshifts of dusty galax- 
ies are significantly worse than they are for normal field 
galaxies which contain less dust. Caution should be exer- 
cised when using large photometric datasets to quantify 
the aggregate properties of infrared-selected samples. 

3.4. Spectroscopic Completeness 

With a lack of any better constraint on the redshift 
distribution of unconfirmed HSGs, we use our targets' 
photometric redshifts (those with and without spectro- 
scopic redshifts) to constrain sp ectr oscopic completeness 
as a function of redshift. Figure 12 (top panel) contrasts 



the spectroscopic and photometric redshift distributions 
by instrument, since LRIS and DEIMOS observations 
have different wavelength coverage, thus different red- 
shift completeness levels. Meauring spectroscopic com- 
pleteness requires an estimate of the underlying redshift 
distribution for the whole population, including those 
without photometric redshifts. The only constraint we 
have on sources without photometric redshifts is from 
those spectroscopically observed, shown in Figure [TT| at 



top, consistent with the distribution in photometric red- 
shifts. Therefore, a scaling correction factor, C, is ap- 
plied to the distribution of photometric redshifts to ac- 
count for the sources without photometric redshifts. Al- 
though this assumption might not be entirely correct 
(i.e. there might be a much higher fraction of z > 2 
sources which do not have photometric redshifts), the 
assumption that the redshift distribution is weighted at 
lower redshifts effectively places the lower limit of spec- 
troscopic completeness at z < 2. If the fraction of sources 
without photometric redshifts is higher at high-z, then 
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Fig. 12. — The top and middle panels show the redshift distri- 
butions of HSGs split into LRIS and DEIMOS samples. DEIMOS 
redshift desert is marked in the middle panel. Assuming the red- 
shift distribution for sources without photometric redshifts resem- 
bles the distribution in photometric redshifts (see Figure |11| top 
panel), then we can estimate the spectroscopic completeness of our 
survey, bottom panel, by dividing the spectroscopic redshift distri- 
butions (by instrument) by the photometric redshift distribution 
scaled up to account for sources without photometric redshifts. 
From this, we determine that survey completeness cannot be con- 
strained at z > 2 (due to small number statistics) which is why 
that sample is discussed in a separate paper. The best-fit curve at 
z < 2 (plotted in gray with dashed lines showing uncertainty) is 
used later to estimate the total contribution of HSGs to the SFRD. 
The two 'outlier' high-completeness DEIMOS points are discussed 
in the text. 



the spectroscopic completeness measurement at z < 2 
will increase. The spectroscopic completeness is mea- 
sured as a function of redshift by dividing the distri- 
bution in spectroscopic redshifts by that of the scaled 
up photometric redshifts (i.e. N(z spcc )/(C x N(z p Y lo t)) )■ 
This spectroscopic completeness estimate is shown in the 
bottom panel of Figure [12] As expected, it declines 
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Fig. 13. — Herschel-SPIKE color with redshift, where color is given 
as a magnitude difference between Spire bands, such that high val- 
ues correspond to red colors. In light gray are Spire colors of the 
entire photometric HSG sample from the COSMOS field; dark gray 
points are our spectroscopic sample. The mean of the spectroscopic 
sample and photometric sample is shown as black circles and dark 
red squares, respectively, with associated uncertainties. Overplot- 
ted are tracks of SEDs with fixed dust temperature, 20 K (red), 
30 K (orange), 40 K (green), 50 K (blue), and 60 K (purple). We 
see no significant color evolution in the HSG sample with redshift; 
this is primarily due to selection properties of the sample. 

over < z < 2 and then is unconstrained at z > 2 
due to limited samples. The DEIMOS completeness es- 
timate has two 'outlier' points at z ~ 0.8 and z ~ 1.3; 
this is caused by the enhanced sensitivity of DEIMOS 
from 6500-7400A (where LRIS often has a gap in cover- 
age) , and 7600-9000A (where our LRIS observations were 
slightly more prone to skyline contamination). Since the 
completeness becomes unconstrained at z > 2, we choose 
z = 2 as a natural boundary at which to divide the sam- 
ple, addressing the well-constrained z < 2 population in 
this paper, while discussing the full sample of spectro- 
scopically incomplete z > 2 sources in C12. 

3.5. Spire Color with Redshift 

Recalling from the discussion in § [31 our spectroscopic 
HSGs sample the underlying distribution of all Spire 
sources well in Spire color, despite a bias in i-band mag- 
nitude (where brighter optical sources are more likely to 
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be spectroscopically confirmed) . Figure [8] showed that 
the mean redshift and Spire color per z-band magni- 
tude of our spectroscopic HSGs is drawn from the parent 
sample of Spire sources without significant bias. While 
this hinted at the evolution of Spire color with red- 
shift, we can test this evolution directly using our spec- 
troscopic HSGs and larger photom etric redshift samples 
from COSMOS flllbert et al. pOlO] ). 
Assuming that most infrared starbursts can be de- 



scribed with a 35 K dust SED (Chapman et al. 2005 



Rieke et al. |2009| ), one would expect low redshift sources 



to have blue Spire colors and high redshift sources to 



2011 Roseboom 



have red Sp ire colors (e.g. |Cox et aT 
et al. |2012[> Figure [l3] plots the [350j-[ - 250j [500j-[350j , 
and |500J-|250] colors against redshift, where [250] is de- 
fined as a 250 /im magnitude, i.e. — 2.51og(S l 2 5o)- We see 
no strong trend in SED shape with spectroscopic red- 
shift. While it might be expected that high-z sources are 
significantly more 'red' than 'blue' low-z sources, this in- 
terpretation is simplistic since it assumes no evolution in 
dust temperature, no relationship between dust temper- 
ature and luminosity, and ignores the impact of popula- 
tion selection effects. Overplotted on Figure [13] are the 
color-redshift tracks for SEDs of fixed dust temperature 
systems from 20-60K; they evolve more strongly with 
redshift than our data, particularly notable in [500]-[250] 
color. 

The dispersion on the Spire colors in both the spec- 
troscopic and photometric samples imply a wide range 
of dust temperatures. This highlights that far-infrared 
flux densities and colors cannot be used exclusively as 
a proxy for redshift. Notwithstanding the large disper- 
sion in Spire colors, two main observations can be drawn 
from Figure 



13 



(a) that the mean Spire colors do not 
redden significantly with redshift as might be expected; 
and (b) at low-z (z < 0.5) our spectroscopic sample is 
redder, i.e. cooler, than the photometric redshift sample. 

Both of these observations can be understood by in- 
vestigating how Spire selec tion works in the Lm — Tdust 
plane, discussed fully in § |4.4| Dust temperature in- 
creases with luminosity (which is partially a selection ef- 
fect against high-z cold systems and partially thought to 
be a real correlation). Since the highest-redshift sources 
in our sample are significantly more luminous than our 
lower redshift sources, the correlation between Ljr and 
Tdust translates to a correlation between z and Td us t- 
This results in roughly constant Spire colors with red- 
shift. At z < 0.5, we note that our spectroscopic survey 
is limited in the number of very luminous ( ^ 1O 12 L0) 
galaxies it can detect simp ly due to sky area probed 
(rj Ideg 2 ); since (in § 4.4) we will observe that more 
luminous HSGs are hotter, it is natural that a sample of 
HSGs with photometric redshifts (covering a much larger 
area than our spectroscopic survey) will probe naturally 
warmer HSGs. This gives rise to the 'redder' z < 0.5 
spectroscopic sample. Although the two redshift bins at 
1.6 < z < 2.0 seem bluer than at any other redshift, we 
note that the number of sources in these redshift bins is 
significantly fewer (and the observation is not the same 
for the much larger sample of sources with photometric 
redshifts in that range). 



4.1. Summary of Spectra and Redshifts 

Table [2] lists the positions, redshifts, Spire flux den- 
sities, 24 /im and radio flux densities, luminosities and 
dust temperatures for the full sample of z < 2 spectro- 
scopically confirmed HSGs. Sample spectra of a variety 
of spectroscopically-confirmed targets of the z < 2 con- 
firmed HSG sample are shown in Figure [14] Of 1594 tar- 
gets with >3cr Spire detections, 767 have spectroscop- 
ically identified redshifts from these data, 731 of which 
are at z < 2. By field, there are 288 in COSMOS, 163 in 
LHN, 119 in GOODS-N, 139 in EN1, 51 in CDFS, and 
8 in the UDS. The total sky area probed by this redshift 
survey is ~0.93deg 2 . 

Besides the selection and spectroscopic biases men- 
tioned in the previous section, the spectroscopic yield for 
sources varied strongly as a function of observing condi- 
tions and integration time, anywhere from 15% to 80%. 
However, unlike the selection and spectroscopic biases, 
this variation in yield is much easier to correct. While the 
full 0.93 deg 2 , 767 source sample is used in most analysis, 
only the sources surveyed under photometric conditions 
to full survey depth, i.e. 13 LRIS masks and 16 DEIMOS 
masks, covering an area of 0.43 deg 2 , are used to measure 
quantities which rely on accurate measurements of source 
density (e.g. the integrated infrared luminosity function 
and star formation rate density). 

4.2. Infrared SED Fitting 

We note that best-fit SEP tem p late fits from|Chary fc| 
| Elbaz|(|200l1),|Dale et aT\ IfgOO lJ) , [Dale fc Helou~i[ 2002 ) , 
|siebenmorgen fc Krugel| ( |2007p or |Draine & Li| ( |2007| ) re- 
cover Lm well, however we choose not to fit our data to 
these complex models given the limitations of our data 
and th e variation in SED shapes. As discussed in |Casey| 
(2012 1 at length, the templates should be used with cau- 



tion when the number of free parameters in the templates 
exceeds the number of data points; this is certainly the 
case for those sources which have three to five photo- 
metric FIR data points. The dust temperature range of 
these SED templates is restricted to the 20< 7d ust <60 K 
range and quantized. For these reasons, we opt to use 
the functional fits in lieu of FIR SED templates. 

Measuring infrared luminosities and dust temperature 
for Spire galaxies is performed by fitting a modified 
blackbody SED extrapolated over the rest-frame FIR 
to existing photometric data. In addition to the Spire 
flux densities, we use photometry from _ffersc/ie/-PACS 
at 100 /im, and 160 /im in Lockman Hole North and 
COSMOS, and Spitzer-MlPS at 24 /mi, and 70 /im where 
available. 

We fit the data for each galaxy to a modified blackbody 
fit of general opacity coadded to a mid-infrared power 
law, where the flux density, S u at rest-frame frequency v 
is represented by 



bb - 
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4. RESULTS 



where t(v)=(u/uo)^ is the optical depth ( Draine|[2006 ), 
j3 is the emissivity, the physical dust temperature is T, 
and v=c/\$ (Aq — 200 /i m) is the fre quency where op- 
tical depth is unity (see Casey 2012 for full details on 
the SED fitting method). Note that the physical dust 
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Fig. 14. — Ten sample spectra of Herschel-SPIRE selected galaxies ranging in redshift from z = 0.126 to z = 1.940. These sample spectra 
span a range of spectral types and include both LRIS and DEIMOS observations. The flux scales are in arbitrary, uncalibrated flux units. 
Spectra of the z > 2 HSGs are given in full in C12. Individual spectra for sources given in Tableware available on request. 

temperature, T, used in this equation is not the same 
as the dust temperature used for the remainder of the 
paper; we choose to measure a "peak-SED" dust tem- 
perature, called Tdust, which is the dust temperature 
measured from the peak of the SED via Wien's law. 

is less dependent on model parameters than T, so is 
more easily comp ared to dust temperatures measured via 
other means (see |Casey|2012l , Fi gure 2 for more detail) . 
Nbh an d -Wpi are the coefficients of the modified black 
body and power law terms, where N p \ is set such that 
the two functions are equal at the frequency, u c , where 
the gradient of the black body is equal to the slope of 
the power law, a. We allow a to vary with the range 
0.5 < a < 5.5 for sources with 24 /jm measurements (all 
sources except those in GOODS-N). If we have no mid- 
infrared measurements for the source, we fix a = 2.0, 
which is the mean value for the rest of the Spire galax- 
ies in our sample and a common value used throughout 
the literature (see Younger et al. 2009 Magnelli et al. 



2010||Hilton et al. ||2012 ). Since our observations cover 
a relatively narrow wavelength range without measure- 
ments on the Rayleigh-Jeans tail in the millimeter (from 
«850 /im-2mm) , we cannot constrain /3 in a meaning- 
ful way (except for very low-redshift sources where the 
Spire points sit on the Rayleigh-Jeans tail). For con- 
sistency in our fitting technique, we decided to fix the 
emissivity to 8 = 1.5 (this is a commonly chosen value of 
B in the literature, e.g. |Chapman et al. |2005{ |Pope et al] 



J|2005J |Younger et al. ||2009j ). There is very little change 
m derived Lir or 'J dust by fixing 8 within the 1-2 range, 
so we determine this procedure to be reasonable. 

Figure [15] illustrates ten example SED fits using our 
method, contrasting them with both radio-implied FIR 
SEDs (from the radio/FIR correlation) and extrapola- 
tions from 24 /im assuming a nomina l template for high-z 



ULIRGs, like the Pope et al. (2008) composite SED 



Table [2] lists the 8-1000 /im 1R luminosities and dust 
temperatures for the whole sample, and Figure fl6] shows 
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Fig. 15. — Ten example FIR SED fits to SPlRE-selected star-bursts: three from LHN; three from COSMOS; two from EN1; and two 
from GOODS-N. The FIR data points (Spi re, and PACS, where available) are shown with ler uncertainties, and our best-fit modified 
blackbody SED is shown in black (see § |4.2| for details on SED fitting). The underlying cold-dust blackbody dominating the FIR portion 
>40 um of our SED is shown a s a dashed line. For sources which have Spitzer-MIPS coverage, we overplot the mean composite SED of 
SMGs from |Pope et al | ( [2008] ) normalised to 24 um (red). Sources with 1.4 GHz VLA coverage have another SEP overplotted (blue): a 
modified blackbody mtered from assuming the radio/FIR correlation holds (with gin, as given inllvison et al. 120101 with T^ ust best fit to 
the data. Gray vertical lines mark rest-frame 8 (mi, 40 (im, and 1000 um. We note that the cold-dust blackbody only dominates the SED in 
the 40—1000 fim wavelength range, although IR luminosity is computed across the range 8—1000 fim, using the combination of IR modified 
blackbody and mid-infrared power law. 

and flux density. Similarly, the radio and 24/im detec- 
tion limits are not absolute, since survey depths at both 
wavelengths vary between fields. 

Figure 16 shows all spectroscopic HSGs in Lir-z space. 
Despite the fact that our sources are identified initially by 
24/im or radio source observations, the detection thresh- 
old in LiR-z space corresponds well with the Spire >3er 
luminosity limit alone (at z < 3, where Spitzer is effec- 
tively more sensitive than Spire). 

4.3. Aggregate Infrared SEDs 

While individual galaxies in our sample have at most 
five flux density measurements in the FIR, some with 
very low signal-to-noise, we can combine the measure- 
ments from many sources to infer the aggregate infrared 
properties of our sample in greater detail. 

Figure [IT] shows two composite near-infrared through 
radio SEDs for the entire < z < 2 HSG sample. The 
first comprises the 331 sources which are radio-detected 
at z < 2: their flux densities are scaled so that the 
i^-corrected radio flux density equals the mean of the 
sample, 133 /iJy, at the mean redshift of the sample, 
(z) = 0.89. The second composite SED is constructed of 
all 731 z < 2 sources, with flux densities re-normalized to 
the mean infrared luminosity, 4.1xlO n L Q , and redshift, 
(z) = 0.80, of the sample. The two SEDs are remark- 
ably similar, with dust temperatures of 27 K and 30 K, 
and mid-infrared slopes of a = 1.9 and 2.1, respectively. 
The most noticable difference lies in the near-infrared, 
where radio normalization seems to artificially wash out 
the stellar mass bump. The other observation to make 
is the spread in the radio flux densities in the IR scaled 
SED; although this is large (~2dex, most falling within 
~1 dex scatter), we note that this could very well be due 




Fig. 16. — The 8-1000/tm integrated infrared luminosity against 
redshift for .fler.se/iei-SPlRE spectroscopically identified sources. 
The lower luminosity limits at the selection wavelengths are 
marked, both for the prior source catalogs at 24 /im (S24 >100/iJy, 
dashed line) and 1.4 GHz (Si. 4 >25 or 40 fi3y, long dashed and 
dot-dashed lines), and Spire wavelengths (colored hashed regions). 
The detection limits at the Spire bands change with SED shape; 
we measure the dust temperatures of galaxies (via Wien's law) se- 
lected in each filter to determine the range of peak-SED dust tem- 
peratures to integrate between to calculate the IR, luminosity limits 
at 250 /^m (blue), 350 fim (yellow), and 500 /im (red). Sources se- 
lected as radio sources are marked as crosses and 24 fim sources are 
marked as green diamonds. 

how the sources fall in luminosity versus redshift with 
respect to the various selection criteria. Note that the 
Spire selection criteria are not fixed in infrared luminos- 
ity since luminosity depends both on dust temperature 
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Fig. 17. — Composite infrared to radio spectral energy distribu- 
tions of HSGs scaled to their radio flux densities (top; 331 radio- 
detected HSGs) and integrated infrared luminosities (bottom; all 
767 HSGs). The best fit SED shapes (light blue, top, dark red, bot- 
tom) are remarkably similar between the two samples, with best-fit 
dust temperatures of 27 K and 30 K and mid-infrared slopes of 1.9 
and 2.1, respectively. 

to individual source variation in the synchrotron slope 
or bright radio AGN. Since the mean value of the radio 
flux densities falls within ±0.2 dex of the expected radio 
flux density (from the FIR radio correlation) , wc deduce 
agreement. The overall agreement between the two ra- 
dio scaled and IR scaled SEDs re-enforces that: (a) the 
FIR radio correlation seems to hold in this sample within 
uncertainty; and that (b) radio selection is not clearly bi- 
ased in integrated infrared properties in comparison to 
24 /im selected HSGs. 

To test for luminosity and redshift evolution of infrared 
SED type, Figure [T8| illustrates the near-infrared through 
radio SEDs of all < z < 2 sources split into several in- 
frared luminosity bins. We measure a steady increase in 
dust temperature with luminosity (and also consequently 
redshift) from «20K at 10 10 L , to «30K at 10 12 L , 
and a jump to 45 K at io 12 - 5-13 L . Understanding 
whether or not this increase in dust temperature is due 
to redshift evolution, luminosity, or selection effects is 
difficult, due to the redshift dependence of luminosity; 
this issue is discussed more in context of the Lir — 7d us t 
plane in the next section. 

An interesting conclusion can be drawn from the 
composite SEDs about the evolution of near- and 
mid- infrared properties. At lower luminosities (up to 

10 12 L ), the near-infrared data <10/im increases to- in this composite SED only due to sources at the corre 
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Fig. 18. — Composite infrared to radio spectral energy distri- 
butions of HSGs in six luminosity bins, from Ljr =10 10 Lq to 
10 13 Lq. Flux densities are renormalized first to the mean lu- 
minosity and redshift of each bin. Infrared dust SEDs are fit to 
binned data (red points) from rest-frame 8 /im to 1000 fim, with 
fixed a = 2.0 and f3 = 1.5. The radio portion of the SEDs are gen- 
erated by assuming the FIR/radio correlation holds; overall, these 
radio SEDs agree with radio data. From the SED fits, we mea- 
sure a steady increase in dust temperature with luminosity, from 
T =19K at 2xl0 10 L Q at (z) = 0.23 to T = 46 K at 5.2xl0 12 L 
at (z) = 1.81. We also see evidence for the evolution of near- to 
mid-infrared properties, most prominent in the highest luminosity 
bin, with detection of the RilO (im Si absorption feature (from a 
24/^m flux density deficit of sources around z ~ 1.4). 



imum at ~5-6 /im. In the highest luminosity bin, the 
prominence of the stellar bump disappears, as the emis- 
sion seems to be dominated by a power law. Further- 
more, the rest-frame mid-IR is unremarkable at low lumi- 
nosities, largely consistent with the expected continuum 
flux density from our SED fit, while the highest lumi- 
nosity bin shows a strong absorption feature at wlO/im, 
thought to be 9.7/iin Si absorption, seen in some of the 
most luminous local infrared galaxies, including Arp 220 
( |Charmandaris et al. 1997[ ). This feature is detected 



wards shorter wavelengths, indicative of emission from 
old stars (i.e. the stellar mass bump) with a local min- 



sponding redshift (z ~ 1.4) having a 24/im-flux density 
deficit relative to sources at redshifts above and below 
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z ~ 1.4. 

4.4. The Temperature-Luminosity Plane 

Figure [19] shows our sample in Tir — Tdust space. Over- 
all, dust temperature increases with luminosity, which is 
partially attributable to our selection and partially to a 
real physical effect. The selection effect (that more lu- 
minous sources are hotter) stems from the nature of se- 
lection at wavelengths on the Rayleigh- Jeans tail, where 
hot-dust SEDs are selected against (at comparable IR lu- 
minosities). This dust-temperature selection effect was 
much more pr onounced in the SMG population selected 
at 850 /mi (see|Blain et al. ||2004HChapman et al. I[2004| 



Casey et al. 20U3~Chapin, Hughes & Aretxaga 2009) 



than it is in Spire samples, which probe the dust SED 
closer to its peak and even beyond, at z ^ 3. The solid 



lines in the left panel of Figure 19 represent Spire de- 
tection limits, or contours of fixed 250 /im flux density 
and redshift. When split into different redshift bins, the 
Lm ~ ?dust relation appears to evolve, although this is 
primarily attributable to luminosity-limit selection ef- 
fects. In order to probe intrinsic redshift evolution of 
sources in Tir— Tdust > a larger dynamic range of luminosi- 
ties is needed over narrow redshift ranges (e.g. expand- 
ing to > 10 12-13 L Q sources at 0.0 < z < 0.5 would help 
greatly) . In contrast to 850 /im selection, 250 /im selec 
tion is relati vely unbiased with dust temperature ( Casey 
et al. |[20TT| ). 



There has been work which argues that the selection bi- 
ases in Ttr — Td ust space ar e minor compared to physical 
effects ( Sajina et al. ||2007 ). The argument is that there 
is more intense dust heating in more extreme star form- 
ing environments. Our data support this, since the dust 
temperature bias effects are less prominant using Spire 
selection. Although selection effects play a role in our 
perceived relation, we fit dust temperature as a function 
of Lir and find a significant correlation for Spire selected 
samples of (T dust ) = 0.47±g;gi (Tfir)°- 144±0 - 006 [K T^ 1 ]. 
The slope of this correlation, 0.14, is steeper than is seen 
in the local IRAS sample. We highlight that the differ- 
ence between the two samples can be attributed to two 
factors. First, the dust temperature bias of the IRAS 
sample, selected at 60 /im with S eo > 5.24 Jy, which se- 



lects against low dust temperature sources (see Figure 19 
right). Second, the local sample is a single-redshift snap 
shot at z = 0, while our sample is spread over a wide 
epoch range, < z < 2 (with 36 sources at 2 < z < 5). 
Fitting the whole sample together washes out any evolu- 
tion in the Tir — Tdust relation d ue to the strong redshif t 



2010) 



dependence of luminosity (e.g. Seymour et al. 
Eventually, much larger samples will make it possible to 
probe the Tir — Tdust relation as a function of redshift. 

4.5. Comparison to SMGs and 'missing' warm-dust 
ULIRGs 

Submillimetre Galaxies (SMGs) selected at 0.85- 
1.3 mm by SCUBA, LABOCA, AzTEC, or MAMBO 
have dominated the studies of high-z starbursts for the 
past decade. While /Jerscfe d-selected galaxies overlap 



with the SMG population ( Roseboom et al. 
the correspondence between source detectability at 250- 
500 /im and ^850 /im is not obvious. Although Herschel 
detects galaxies at their SED peak at z ~ 1 — 2 (and 
is more sensitive to lower luminosities at z < 2 than 



SCUBA), the submm has the advantage of preferentially 
detecting higher-redshift galaxies and detecting them to 
lower luminosity limits than Herschel. 

Unfortunately, very few of our Spire galaxies were 
directly observed with SCUBA, LABOCA, AzTEC or 
MAMBO on account of small survey areas for those in- 
struments, so it is difficult to assess population overlap 
directly. 



4.2 



However, using the SED fits from 
can estimate 850 /im flux densities for every HSG. These 
SED-predicted 850 /im flux densities span a wide range, 
from ^0.2-30 mJy. Figure 20 shows the relation between 
fitted dust temperature and extrapolated 850 /im flux 
density for Spire galaxies; there is a clear relation be- 
tween dust temperature and 850 /im-detectability. 

We infer that 79lg% of all Spire galaxies are unde- 
tectable at 850 /im at a flux cut-off of Ssso < 5 mJy. Even 
considering a submillimeter detection threshold as low as 
2 mJy (the lowest 3tr detection limit for SC UBA, which 
had a confusion noise of ^0.7 mJy at 850 /im Blai n et al.\ 
1999), 47+H% 01 Spire galaxies would be undetectable 
at 8t>0 /im. While this is a large fraction, it could be 
that most of these are at low redshift, and therefore low 
luminosities and low SFRs. The 'missing' fraction as a 
function of redshift is plotted on the right panel of Fig- 
ure [20j Indeed the fraction of HSGs not detectable at 
850/mr is high at z < 1, 83+7%, yet even at higher red- 
shifts, 1 < z < 2, the fraction is substantial: 727_i!% of 
HSGs have Ss5o < 5 mJy. 

The fact that 850 /im selection misses a large fraction 
of infrared-bright starbursts is not new, but has been dif- 
ficult to measure directly or estimate in the past. This 
dust temperature selecti on effect was studied in detail 
in the pre-Herschel era dBlain et 



al. 



eFar|[2004l [Casey et al. | |2009[ |20ll|? |) 

tive estimates of the 'missing' warm-aus 



2004 
whe 



Chapman 



ire conserva- 

tive estimates of the 'missing' warm-dust ULIRG popula- 
tion (often refered to as 'OFRGs' or Submillimeter-Faint 
Radio Galaxies, 'SFRGs') were on the order of 10-20%, 
and later works went on to estimate the missing fraction 
at £ 30% using photome tric data and limited spectro- 
scopic data from Herschel (Chapman et al. 2010 Magdis 



et al. 



2010). 



While we note that roughly half of Spire galaxies are 
below the detection threshold at 850 /im, we can also in- 
vestigate the converse statistic. How many SMGs are 
not detectable by Spire? We use the ori ginal sample of 
SMGs described in |Chapman et a l. ( 2005 ) , along with as- 
sumed IR luminosities estimated via the FIR/radio cor- 
relation to fit SEDs and extrapolate to Spire fluxes. We 
then use these hypothetical fluxes to calculate a Spire 
detectability, in other words, we determine what fraction 
of the SMGs would pass our >3er selection criteria. The 
percentage of SMGs which are 250-500 /im detectable is 
shown at bottom on the right panel of Figure [20] At 
z < 2, all SMGs are detected with Spire, but above 
z > 2, anywhere from 20% to 60% of SMGs are Spire 
dropouts. 

We also investigate this fraction using real Spire cov- 
er age of SMG fields togeth er with LHN data presented 

and Spire coverage around 
'( [20051 ). 



Chapma n et al. (2010) 
2012D , SMGs irom|Chapman et al. 



The statistic is con- 
sistent, roughly 1/3 of SMGs are not detected in Spire, 
and nearly all at z > 2.5. Since the sample size of SMGs 
with overlapping Spire data is about half the size of the 
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Fig. 19. — Lir against 7d uat as an indicator of the variation of SED type in our sample. Dust temperature is estimated via Wien's Law 
(inversely proportional to peak of the SED). Colors on the left denote redshift bins, and the mean Ljji — Td us t relation for each rcdshift 
slice is shown as large squares. This relation changes with increasing redshift, however this can be attributable to the selection bias of 
Spire selection: the 250 fim detection limits at fixed redshifts are shown as solid lines, whose color corresponds to the lower redshift limit 
of the bin. At a fixed temperature, only sources to the right of the line (i.e. at higher luminosities) are detectable. The mean luminosities 
and temperatures for the composite SEDs in Figure |T8] are large crosses. At right, we illust rate the Ltr — T^„„t. relation for Spire ga laxies 
(black) in co mparison to local IRAS galaxies, from the RBGS and GOALS samples (red; [Sanders et al. |2003| |Armus et al. ||2009| l, and 
SMGs (blue; Chapman et al. 2005; Pope et al. ||2006] l. We refit SEDs for all samples consistently so that both dust temperatures and 
luminosities are directly comparable. This demonstrates that our Spire galaxy sample is statistically colder than local IRAS galaxies of 
similar luminosities but warmer than SMGs of similar luminosities. Although we observe a different slope to the Ltr — Th ,,^ relation than 
is seen locally, our sample is not large enough to measure evolution, since luminosity is largely a function of redshift (e.g. Seymour et al. \ 
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Fig. 20. — Predicted 850 fim flux densities for our Spire sample against the fitted SED dust temperatures. The characteristic uncertainty 
in either measurement is shown in the upper right. Note that the classic SMG selection criteria, Ssso > 5mjy would exclude 79^g% of 
Spire sources from detection in traditional 850/xm-based SMG surveys: 83^% of HSGs at z < 1 and 72^% of HSGs at 1< z < 2. This 
is significantly higher than previous estimates to the 'missing' fraction of warm-dust ULIRGs at z > 1 missed by SMG surveys. At right, 
we show the fraction of Spire galaxies which are 850/xm-detectable and 850^tm-selected SMGs which are SPlRE-detectable, as a function of 
redshift. Here Spire detectability is defined as S25O1 or 6*500 > 12mjy, and qualification as an SMG is Sgso > 5mjy. This illustrates 

that 60-80% of Spire galaxies are undetectable at 850 fim out to 2 ~ 1.6, and that SMGs are only expected to drop out of the Spire bands 
at 2 > 2 for 20-60% of sources. The numbers next to each point represents the number of galaxies used in each bin. There are too few 
Spire galaxies to reliably determine this statistic at z > 1.6 (gray shaded region). 
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Fig. 21. — Ratio between FIR luminosity and radio luminosity in 
our sample, grR, against r edshift. We compare against the moder- 
ate evolution measured by Ivison et al. |2010[ | oc (1 + z)~ 0- 26 ± 007 
(gray band). We measure slightly stronger evolution, qm oc 
(l + z)-°- 30 . 

full SMG sample, we use the extrapolated flux densities 
in Figure [20[ 

4.6. The FIR/Radio Correlation 

Since IR luminosities were derived independently of 
radio luminosity, we measure the FIR/radio ratio, qm, 
for radio-detected HSGs in ord er to assess the FIR/radio 
correlation for starbursts (see 



Robinson 
'i uh/i _ r 


1985 


Condon 


1992 



Hclou, Soife r fc Row an 
betweer 



vison et al. 



for 

jet ween 8- 



1000 /xm flux and 1.4 GHz radio flux as defined in Ivison 



et al. 



(2010), with IR luminosities measured from rest 
frame 8/iim to 1000/im. Figure [21] shows the qir, ratio 
evolving with redshift for synchrotron slope a — 0.75, 
where S v oc v~ a . Our sample is consistent with the 
mo derate evolution measu red by |Ivison et al. | ( 2010|? 



and Magnelli et al. (2010). We measure qir oc (1 + z) 
where 7 =-0.30±0.02HIsmg only the < z < 2 HSG 
sample (an indepe ndent assessm e nt of the z > 2 sample 
is given in C12). Ivison et al. (2010) measured 7 =- 
0.04±0.03 at z = 0-2, but found thaTtheir z < 0.5 sam- 
ples contaminated this measurement since they were few 
in number and not well matched in luminosity to the 
higher redshift sources; they measure 7 =-0.26±0.07 for 
the 0.5 < z < 2.0 sample, which is consistent with our 
finding, 7 =-0.30. Either with or without our z < 0.5 
sample, we measure 7 =-0.30, likely because our sam- 
ple is dominated (68%) by sources at 0.5 < z < 1.5. 
Note that the evolution we measure is much more pro- 
nounced if the original luminosity limits defining Sfir 
are used, 40-120^m: q FIR oc (1 + z )-o.59±o.04^ Note also 
that the significantly limiting factor of this measurement 
is that our sample is partially radio-selected and is not a 
luminosity-matched sample. 

5. DISCUSSION 

The characterization of galaxies by Herschel is the 
largest advance in understanding star formation in the 
infrared since the discovery of SMGs by SCUBA. This 



paper has presented a spectroscopic survey of 1594 Spire 
sources, 767 have identified spectroscopic redshifts from 
< z < 5, and 731 at z < 2. Securing spectroscopic 
redshifts is itself very valuable for follow-up studies of 
this population, from their metalicities, stellar popula- 
tions, morphologies and evolutionary histories to their 
molecular gas properties and dust content. However, 
arguably the impact of this population on the cosmic 
star formation rate density (SFRD) is the most relevant 
computation for galaxy formation and evolution studies 
in general. The luminosity function has only ever been 
measured in integrated IR luminosity (Ltr(8— 1000/wi)) 
for li mited populations of SMGs with known inco mplete- 
ness ( Chapman et al. 2005 Wardlow et al. |2011 ) or from 
extrapolations to the IR from Svitzer jLe Floe h et al. 



2005| |Caputi et aT"||2007[ |Magnelli et al. ||201ip inak" 



mg use of model templ ates or explicitl y assuming SEDs 



from libraries ( Le Borgne et al 



2011 Marsden et 



2009 Bethermin et al. 



2011). This section presents the 



integrated IR luminosity function of Herschel galaxies, a 
discussion of the completeness of our samples, the impli- 
cations for the cosmic SFRD and comparisons to these 
previous studies. 

5.1. Luminosity Function 

We compute the integrated IR (8-1000 ^im) luminos- 
ity functions (LFs) of SPiRE-detected galaxies in four 
redshift bins where the spectroscopic completeness is 
constrained and the number of identifications is ^100 
per bin: 0.0< z <0.4, 0.4< z <0.8, 0.8< z <1.2, and 
1.2< z <1.6. We use the 1/V max accessible volume tech- 
nique where 

™ 1 

$ 2 (i)AL = V^^ (2) 

$ Z (L) is given in units of h 3 Mpc -3 logL -1 and $ Z (L) AL 
is the number density of sources with luminosities be- 
tween L and L + AL. Here Vl(L) is the comoving vol- 
ume within which the zth source is detectable in our sur- 
vey. This accessible volume is determined by construct- 
ing a detection limit in Lir—z space for each source. As 
was evident in Figure 16 the Spire detection limit is 
highly dependent on dust temperature (i.e. 40 K galax- 
ies will have a different -Ltr-z detection threshold than 
20 K galaxies). The luminosity limit per source is de- 
termined by assuming an SED (of type given by Eq 1) 
with fitted Tdust , and the 3<r detection bounds at 250 /im, 
350 /ini, and 500 /im, where the lowest luminosity limit 
with redshift is adopted for the accessible volume calcula- 
tion. The measured luminosity function, after correction 
for spectroscopic incomplete ness (as discussed earlier in 
the text and shown in Figure 12 ), is given in Table [l] and 
shown in Figure 22 as filled, colored circles. 

As a check, we also compute the luminosity func- 
tions for the same XID se lection using the p hotomet- 



2010| ) and the 
^6000 sources over 



ric redshifts in COSMOS ( jllbert et al. 
same luminosity limit technique (for 
^2.2 deg 2 ). The photometric redshift results are shown 
on Figure [22] as open gray squares. Within uncertainties, 
the spectroscopic and photometric luminosity functions 
agree at each epoch. 

Figure [22] clearly demonstrates evolution in the lu- 
minosity function with redshift, with ultraluminous 
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(>10 L Q ) and hyperluminous (>10 L©) galaxies be- 
coming far more abundant with increasing z than in 
the local Universe. At each of the four epochs, we fit 



Schechter functions of the form $ = 



where 



Fig. 22. — Evolving luminosity functions for our HSG spectroscopic sample. The redshift bins are 0.0 < z < 0.4 (purple), 0.4 < z < 0.8 
(red), 0.8 < z < 1.2 (yellow) an d 1.2 < z < 1.6 (green). We also include the estimated luminosity functions from photometric redshifts 
in COSMOS filbert et al. |2010[ ) as gray boxes and associated uncertainties. The best-fit Schechter functions are shown as dotted colored 
lines, with associated uncertainties (shaded regions) based on uncertainty in a, <!>* and L* . The lower luminosity completeness limits of 
each bin are shown as dashed vertical lines. Below these luminosities, our Spire samples are incomplete ( in t erms of Spire detectability). 
Note that these luminosity functions are corrected for spectroscopic incompleteness (as detailed in Figure [T2"} . 

frared luminosity; while a double power law luminosity 
function might provide an alternate fitting method, the 
difference in this work is inconsequential considering the 
uncertainty on the luminosity function points themselves. 

The parameters $* and L* are then varied to provide 
an optimal fit to data at each epoch. We measure evolu- 
tion in both parameters, as given in Table [l] However, 
it should be stressed that the measurement of these pa- 
rameters' values is highly degenerate with variation in 
the faint end slope, a and with each other. Much larger 
data sets, particularly those which constrain the faint 
end slope via stacking techniques, are needed for detailed 
analysis on the evolution of <&* and L* . This detailed 
work is being done with larger samples from COSMOS 
(M. Vaccari, E. Le Floc'h, private communication). 



x = L/L*. The faint end of the luminosity function, 
dominated by the x a power law, is largely unprobed by 
our dataset since L* , the turnover luminosity, is approx- 
imately equal to the lower luminosity limit of Spire in 
each bin. Without constraints, we decided to fix a =- 
0.35, our measurement from the local RBGS sample 



(Sanders et al. 2003) and the < z < 0.4 photometric 



redshift sample, which is a lso consistent with prior local 



and high-z m easurements ( Le Floc'h et al. |2005 Caputi 



et al. 



2007 ). We add a large margin of uncertlanty on 
the slope a, since its impact on integrated star forma- 
tion density increases with redshift (a a =0.1 at z 0, 
increasing to 0.5 at z rs 1.4). At the top end of the lumi- 
nosity function, we note slight excesses at 0.0 < z < 0.4, 
0.4 < z < 0.8 and 0.8 < z < 1.2, which, in other sam- 
ples, has been attributed to AGN contributions to in- 



5.2. Star Formation Rate Density 

To estimate the contribution of SPiRE-detected sources 
to the cosmic star formation rate density (SFRD), our 
measured luminosty functions are converted to SFR func- 
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TABLE 1 

IR Luminosity Function from the Herschel-SpmE 

SPECTROSCOPIC SAMPLE 
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of galaxy in this context, it is important to point out 
that this depends entirely on the adopted definition of 
'SMG' and associated luminosity limits. The ~ldex 
difference between the SFRD contributions of 850/im- 
selected SMGs and 250-500 Lsm-selected HSGs at z ~ 1 is 
due to the effective luminosity limits of Spire or 850/im 
selection (where HSG selection probes luminosities ~5 
times fainter) and also due to the dust temperature bias 
of 850/im-selection (another factor of ^2). 

Drawing a fair comparison with other samples requires 
interpolation in the luminosity function across equal lu- 
minosity bins via our best-fit Schechter functions. Sev- 
eral literature sources have estimated the total integrated 
star formation rate density from infrared sources using 
a variety o f datasets: Spitzer-24 / xm-identified sources in 
COSMOS flLe Floc'h et al. |2005|); Spitzer-24 /xm sources 



in both GOODS fields (|Caputi et al. | |2007jjMagnelli 



et al. 



2011); and ^fcari-seiected sources (Goto et al. 



We also dra w on the full optical and u ltraviolet 



data compilation of Hopkins & Beacom (20061, which is 
extinction corrected to account for infrared emission in 



Note. — The luminosities in this table are integrated within 
8-1000 /im. They represent the central luminosity of sources in 
each bin, in other words, the bin for luminosity L covers sources 
with luminosities between L — AL/2 and L + AL/2. The lumi- 
nosity functions in this table are based exclusively on our Spire 
galaxy spectroscopic sample. We compute the luminosity func- 
tion in units of h 3 Mpc -3 logL - . 



optically-bright galaxies. On the right panel of Figure 23 
we show the breakdown of the SFRD in 1 dex lumi- 
nosity bins, <10 n L o (sub-LIRGs), 10 11 < z <1O 12 L 
(LIRGs), 10 12 < z <1O 13 L (ULIRGs), and total inte- 
grated infrared. On a whole, we observe the same trends 
as other literature works: the total integrated infrared 
SFRD contribution is comparable to the total optical 
and ultraviolet contribution (i.e. just as much star for- 
mation is obscured as is unobscured) LIRGs dominate 
at z « 1, and ULIRGs become increasingly important 
at z > 1. However, a few subtle differences stand out. 



tions via the Kennicutt ( 1998 1 prescription assuming a 
non-evolving Salpeter IMF. JNote that there is some re- 
cent discussion, particularly in relation to very high lu- 
minosity sources at z > 2, that the Kennicutt conversion 
from infrared l uminosity to SFR migh t not hold due to 



For example , both Le Floc'h et al. (2005) and Caputi 
et al. (2007) find higher contrib utions from ULIRGs at 
z < 0.5 than is measured in the Magnelli et al. (20111 
sample or in our sample. Since both prior samples were 
24/im-selected, this could be indicative of the disagree- 
ment betw een 24/xm and total int regrated luminosity. At 
z > 1, the Magnelli et al. (2011) sample shows a deficit 
in ULIRG s with respect to our sam ple (also seen in the 



models of Bethermin et al. 2011). While this deficit 



paper we use it for consistency with previous work on 



the infrared SFRD (Le Floc'h et al. 


2005; Hopkins & 


Beacom|2006 Caput! et al. 120071 IMag 


nelli etal. 201l|. 



ure [23} The left panel shows the raw conversion of the 
luminosity function data points to SFRD contribution, 
including all sources above and below the completeness 



IMF variation (|Swinbank et al. ||2008|); however, in this could be due to SED assumptions, it could also be due 

to cosmic variance and the size of the GOODS fields used 
in Magnelli etal. (e.g., as is known t o be a problem in 
CDFS, which has a submm deficit, see |Weifi et al. |2009| ). 

The uncertainties on the SFRD measurements are 
dominated by the uncertainty in the faint-end slope of 
the luminosity function, which is not well constrained. 
At higher redshifts L* increases; this impact of the faint- 
end slope uncertainty then grows, since sub-LIRGs and 
LIRGs arc not constrained with our data beyond z >1.2. 
This is not only true for our data set, but any data set 
(selected via 24 /im or radio with a PACS or Spire de- 
tection) which does not employ stacking analysis to con- 
strain faint sources. Further work is needed on stacking 
24 /im- and radio-selected samples of galaxies to probe 
fainter luminosities out to high redshifts over large areas 
of sky less prone to cosmic variance. 

It is clear from this work that the importance of 
infrared-luminous star formation in the context of cos- 
mic star formation is high, and equally high as rest- 
frame UV/optical estimates, even though prior surveys 
of SMGs seem to imply that infrared starbursts are rare 
even at z ~ 1-2. Herschel-SPIKE has enabled us to di- 



luminosity limit (dashed lines on Figure 22 ), and without 
extrapolation to lower or higher luminosities. Although 
we excluded 1.6 < z < 2.0 sources from the luminos- 
ity function calculation due to small numbers, we bring 
them back for this computation. The measured HSG 
SFRD appears to increase from z = to z = 1, and then 
fall sharply at z = 1.5; however, this is due entirely to the 
luminosity limits of the sample at the given epochs. For 
example, the highest redshift bin, z w 1.8, has a lower lu- 
minosity limit of logLiR = 12.3, significantly higher than 
the lower redshift bins. Similar luminosity limit restric- 
tions exist for the literature values of the SFRD found 



from SMGs in |Chapman etai\ ( 2005 1 , Wall, Pope fc 
Scottl ( [2008l ) and lWardlow et al. [(jzOTTj ). 
Although SMGs are thought of as a very rare class 
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Fig. 23. — Contribution of Spire galaxies to the cosmic star formation rate density from z = to z = 2. We compare against 
the compilation in Hopkins & Bcacom 2006 of optical/UV SFRD estimates which are dust-corre cted (large grey a r ea, re prese nting the 
associated un certainty of UV/optical estimates). We also overplot the SFRD points estimated by |Chapman et al. | ( [2005} and |Wardlow| 
\et al. | ( |20 1 1 [ ) for SMGs which are substantially lower by ~10x than the optical/UV dust-corrected SEKD and riSGs at z <1.5. At lert, 
we measure the SFRD from this spectroscopic sample by converting the m eas ured luminosity functions to star formation rate functions, 
then integrating. At right, we use the best-fit Schechter functions in Figure [22] to extrapolate over luminosities not directly probed by our 
survey, to measure the total infrared (green), sub-LIRG (blue), LIRG (orange) and ULIRG (red) contributions to the SFRD (colors are 
con sistent across diff erent data sets). The right y-axes give the value of the IR luminosity density which translates direc tly to SFRD via 
the |Kennicutt| ( [199 81 scaling. We gather simi l ar com putations from the literature for comparison: Lc Floc'h et al. ( 2005 I Spitzer samples 
in GOSMUy (hashed regions), |Caputi et al. | | |2007| l Spitzer samples in GOODS (asterisks), and Magncll i et al. | ( |20lip Spitzer samples in 
GOODS (diamonds). 



rectly constrain the far-infrared unlike previous work in 
the area, extrapolating from 24^m. Our observations 
confirm the importance of dust obscured star formation 
across the first 10 billion years of the Universe's history. 

6. CONCLUSIONS 

This paper has presented a large spectroscopic survey 
of galaxies selected in the Herschel-SPIRE 250-500 /im 
bands and followed up spectroscopically on Keck LRIS 
and DEIMOS. Out of 1594 spectroscopic targets span- 
ning 0.93 deg 2 in multiple deep legacy fields, 767 sources 
had identifiable redshifts from their rest-frame UV or op- 
tical spectra, 731 of which are at z < 2. We present the 
following conclusions: 

• The redshift distribution of SPiRE-selected 
spectroscopically-confirmed galaxies peaks at 
z = 0.85 and has a tail of sources extending out 
to z « 5. The vast majority of this spectroscopic 
sample (731/767« 96%) is at z < 2, mostly 
spectrally identified by [Oil], [OIII], H/3, or Ha 
emission. 

• Our spectroscopic sample excludes sources without 
24 /jm or radio 1.4 GHz counterparts. A negligible 
fraction of SPiRE-sources in deep legacy fields are 
thought to drop out in the mid-IR and radio at 
z < 2, meaning our Spire targeting is close to 
complete at z 2. 

• Only ^50% of spectroscopic targets yield redshifts 
(767 with identifications, 826 without). Under pho- 
tometric observing conditions, ^60% yield identifi- 
cations. We constrain the spectroscopic complete- 
ness at z < 2 to 20-80%, depending on redshift, 
using large catalogs of photometric redshifts. The 



sources without spectroscopic identifications either 
were observed in poor weather conditions, have no 
optical counterpart (i.e. no source picked up on the 
slit after 1-2 hrs integration) , or have spectra not 
easily identified with a single redshift (e.g. contin- 
uum without emission or absorption features). 

• We measure the accuracy of HSG 
optical/near-infrared photometric redshifts at 
Az/(1 + z spcc )=0.29, a factor of ^3-4 times 
worse than photometric redshifts for normal field 
galaxies. We determine that the lack of reliability 
in photometric redshifts for HSGs is due to more 
significant optical obscuration and higher overall 
star formation rates (thus higher line-to-continuum 
ratios). We caution that future analysis of aggre- 
gate properties of infrared-selected samples like 
HSGs might be biased by only making use of large 
photometric redshift catalogs. 

• We observe a correlation between dust temperature 
and infrared luminosity which is partially a selec- 
tion effect and partially thought to be physically 
real. The selection effect is caused by the depen- 
dence of Spire flux density on dust temperature 
and the physical motivation comes from the local 
observation that more luminous galaxies have more 
compact, clumpy and therefore hot dust. 

• Spire color does not evolve with redshift in our 
spectroscopic sample, due to the redshift depen- 
dence of luminosity and the luminosity dependence 
on dust temperature (a combination of selection ef- 
fects and physical mechanisms). 

• We infer the aggregate properties of HSGs by com- 
bining the near-infrared to radio measurements by 
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selection method and luminosity. We find little 
difference between the bulk infrared properties of 
radio-selected galaxies and 24 /mi-selected galax- 
ies, and good overall agreement with the FIR/radio 
correlation. Consistent with our results from indi- 
vidual sources, we see an increase in dust temper- 
ature with luminosity in the composite IR SEDs. 
At luminosities >1O 125 L0, we detect the «10/im 
Si absorption feature in the composite, resembling 
the SED of the local ULIRG Arp 220. 

• By extrapolating our infrared SED fits to 850/im, 
we determine that only 79^ 8 % of HSGs would be 
detectable as 'classic' SMGs (Sgso > 5mJy). This 
fraction is highest at low redshift, 83^% at z < 1, 
and still significant, 72+^%, from 1 < z < 2. 

• For the 331 radio-selected sources with spectro- 
scopic confirmations, we measure the evolution in 
the FIR/radio correlation and find qir = (1 + 
2; ^-o.30±o.02 Thjg evolution is stronger than some 
previous measurements, however the sample is bi- 
ased by being radio-detected, thus the intrinsic evo- 
lution of giR with redshift is likely to be shallower. 

• We compute the luminosity functions for HSGs 
in four redshift bins across < z < 1.6 and 
find agreement with predictions from photomet- 
ric redshifts. The luminosity functions are well- 
fit by Schechter functions, with evolving pa- 
rameters (3.2xl0 3 fr 3 Mpc~ 3 at z = 0.2 to 
3.1xl0- 4 fr 3 Mpc- 3 at z = 1.4) and L* (5xl0 10 L© 
at z = 0.2 to 2x 10 12 L Q at z = 1.4). The faint-end 
of the luminosity function is not well probed by our 
data and would only be accessible from Herschel 
data through stacking analyses. 

• We find the HSG contribution to the cosmic 
star formation rate density to be substantial at 
z < 2, essentially comparable to the star for- 
mation rate density from optical/near-IR surveys 

> lxl0~ 2 M Q yr^ 1 h 3 Mpc" 3 . Without relying on 
extrapolations from the mid-infrared or model tem- 
plate libraries, we measure the contribution from 
LIRGs (10 11 < L <10 12 L Q ) to the SFRD peaks 
at z ~ 1, and the ULIRG (>1O 12 L ) contribution 
increases with z, nearly surpassing the LIRG con- 
tribution at z ~ 1.4. 

This work has demonstrated that infrared-luminous 
galaxies form an integral part of galaxy formation across 
a wide range of epochs, particularly at < z < 2. With 
spectroscopic redshifts in hand, follow-up studies can be 



more efficiently carried out to determine the physical and 
evolutionary origins of HSGs as a function of infrared 
luminosity, dust temperature, dust mass, and AGN con- 
tent. This work will then begin to shed light on the dom- 
inant mechanisms producing obscured star formation in 
the Universe. 
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